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Abstract 

In recent years, asymptotic approximation schemes have been de¬ 
veloped to describe the motion of a small compact object through a vac¬ 
uum background to any order in perturbation theory. The schemes are 
based on rigorous methods of matched asymptotic expansions, which 
account for the object’s finite size, require no “regularization” of diver¬ 
gent quantities, and are valid for strong fields and relativistic speeds. 
Up to couplings of the object’s multipole moments to the external 
background curvature, these schemes have established that at least 
through second order in perturbation theory, the object’s motion sat¬ 
isfies a generalized equivalence principle: it moves on a geodesic of 
a certain smooth metric satisfying the vacuum Einstein equation. I 
describe the foundations of this result, particularly focusing on the 
fundamental notion of how a small object’s motion is represented in 
perturbation theory. The three common representations of perturbed 
motion are (i) the “Gralla-Wald” description in terms of small devia¬ 
tions from a reference geodesic, (ii) the “self-consistent” description in 
terms of a worldline that obeys a self-accelerated equation of motion, 
and (hi) the “osculating geodesics” description, which utilizes both (i) 
and (ii). Because of the coordinate freedom in general relativity, any 
coordinate description of motion in perturbation theory is intimately 
related to the theory’s gauge freedom. I describe asymptotic solutions 
of the Einstein equations adapted to each of the three representations 
of motion, and I discuss the gauge freedom associated with each. I 
conclude with a discussion of how gauge freedom must be refined in 
the context of long-term dynamics. 
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1 Preamble and survey 

Consider a small object moving through a curved spacetime. What path 
does it follow? At the level of undergraduate physics, the answer is satisfy- 
ingly simple: if the object is sufficiently small and light, it can be idealized 
as a test particle, which does not affect the geometry around it and which 
moves on a geodesic of that geometry. But what if we do away with that 
idealization? The real object does perturb the geometry around it; how, 
then, does the object move in that geometry, which it itself affects? 

In Newtonian gravity, we may ask the analogous question—how does 
a massive body move when it contributes to the gravitational field around 
it?—and we may happily provide an answer without yet leaving undergrad¬ 
uate physics: if the body is sufficiently spherical, it can be treated as a point 
particle located at its center of mass, and the motion of that center of mass 
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is governed by the external gravitational fields produced by all other masses 
in the system; the object does not feel its own held, at least in so far as its 
center-of-mass motion is concerned. 

However, in general relativity, the situation becomes radically more com¬ 
plicated. Because of the nonlinearity of the Einstein equations, an extended 
object cannot, in general, be modelled as a point particle without invoking 
post-hoc regularization procedures [HE]. Furthermore, in a curved back¬ 
ground, an object, even an asymptotically small one, does feel its own held, 
for reasons discussed below (and elsewhere in these proceedings). Hence, 
the held is said to exert a gravitational self-force on the object. 

One might wonder if this effect is actually relevant in any realistic sce¬ 
narios. The answer, unequivocally, is “Yes”. Suppose we are interested in 
a bound binary of widely separated compact objects of comparable masses 
mi and m 2 moving slowly in each other’s weak mutual gravity. Each of the 
objects is small compared to the other scales in the system (for example, the 
typical orbital separation R, or the radius of curvature of mi’s held at m 2 ’s 
position). In a Newtonian approximation, mi is subject only to m 2 ’s New¬ 
tonian held, which at the position of mi exerts a force (per unit mass) Fj\f ~ 
m 2 /R? ■ But if one requires anything more accurate than the Newtonian ap¬ 
proximation, the moment one steps to the post-Newtonian level, self-forces 
can no longer be ignored: mi’s held, which we can think of as scaling like 
mi/R, modihes the Newtonian helds, giving rise to a post-Newtonian 
force per unit mass that scales like Tpiv rsj ^nFn rsj mim 2 /R? [3|, similar 
in magnitude to any other post-Newtonian effect 0 

In these systems just described, however, each object, while small com¬ 
pared to the radius of curvature of the ambient external held it hnds itself 
in, is not small compared to the other object. What if we are interested in a 
case where there are truly only two scales? Take a binary system of compact 
objects of mass m and M satisfying m <C M, and specihcally focus on the 
regime in which the orbital separation R is of order M. In this case, it seems 
the gravity of m must certainly have a very small effect on its own motion, 
and it must very nearly follow a geodesic of the metric of M. And yet, even 
in this scenario, m’s gravitational self-force cannot be neglected. Although 

^In a certain sense, an object’s mass affects its own motion even in a Newtonian binary. 
Each object follows a Keplerian orbit about the system’s center of mass, not about the 
center of the other object. Since the center of mass is shifted by the object’s own mass, the 
object affects its own motion; more plainly, mi influences its own motion by influencing 
that of m 2 . This is a more indirect effect than the type described above, but in practice, 
distinguishing it from any other post-test-body effect is nontrivial. See Ref. [4] for a 
discussion. 
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the effect is very small over a few orbits, the system continually radiates en¬ 
ergy in the form of gravitational waves (or equivalently, the self-force does 
negative work), causing the orbit to shrink, and eventually causing m to 
collide with M (or plunge into M, if M is a black hole). In other words, 
the self-force has long-term, secular effects on the motion which make it 
impossible to ignore. 

Both of these two types of systems—binaries of widely separated bodies 
whose mutual gravity is weak, and binaries of objects with very dissimi¬ 
lar masses, called extreme-mass-ratio inspirals (EMRIs)—are of increasing 
relevance to modern astrophysics. The prospect of directly detecting grav¬ 
itational waves emitted from compact binaries, and extracting information 
about the binaries’ strong-held dynamics from those waves, has spurred an 
international effort to study them. In the case of widely separated bodies 
of comparable mass, the main method of analysis has been post-Newtonian 
theoryj^ a historied subject modern overviews of which can be found in the 
review articles 13 E] and the recent textbook [7|. In the case of EMRIs, 
the main method of analysis has been self-force theory; for summaries of 
efforts to model EMRIs, I refer the reader to the reviews pm, the more 
up-to-date but brief survey ini, and the contributions of Babak et al. m 
and Warded [l2] elsewhere in this book. 

In this paper, I seek to provide a single, unihed description of gravi¬ 
tational self-force theory. Roughly speaking, this formalism consists of a 
perturbative expansion in powers of the small object’s mass m. Although I 
will refer to EMRIs to motivate many of the methods and problems I dis¬ 
cuss, my focus will instead be on foundational issues related to the problem 
of motion of a small object. My aim is to complement extant reviews Pi 
by concentrating on three themes given limited attention in those reviews: 
(i) self-force theory at arbitrary perturbative order, (ii) differing ways to 
represent perturbed motion, and the asymptotic expansions of the metric 
corresponding to each, and (iii) the relationship between perturbed motion 
and gauge freedom. My presentation mostly follows the methods and view¬ 
point of Refs. [IMS], though it takes additional inspiration from the work 
of Gralla and Wald [2DH12I and the classic papers of Mino, Sasaki, and 
Tanaka |23] and Detweiler and Whiting |24[ I25j . To avoid excessive length, 
rather than striving for complete self-containment, I will often refer the 
reader to Refs. pnadziiiHi for mathematical tools and technical details. 

^Of course, once the two bodies are sufficiently close to each other, they interact in a 
highly nonlinear, highly relativistic way. In that regime, one must use numerical relativity 
to solve the fully nonlinear Einstein equations for the system. 
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Most of what I discuss could be applied to objects carrying scalar or 
electric charges, but for simplicity I restrict my attention to the purely 
gravitational case. I refer the reader to Refs. [9l [26H30] for discussions of 
self-forces due to scalar and electromagnetic fields on fixed background ge¬ 
ometries, and to Refs. [m [32] for recent work on the coupled system in 
which the metric reacts to both the mass and the scalar or electromagnetic 
field. 

In the remainder of this introduction, I make an extensive survey of the 
main concepts and results of self-force theory, beginning with the standard 
picture of a point mass in linearized perturbation theory and proceeding to 
describe the more robust framework now available for studying the motion 
of extended (but small) objects at any order in perturbation theory. 

Notation. Throughout this paper, Greek indices run from 0 to 4 and are 
raised and lowered with a background metric g^i,. Latin indices run from 1 
to 3 and are freely raised and lowered with the Euclidean metric 5ij. Sans- 
serif symbols such as refer to tensors on the perturbed spacetime rather 
than on the background. A semicolon and V refer to a covariant derivative 
compatible with g^,^, to the covariant derivative compatible with g^,^, 
and to the covariant derivative compatible with an “effective metric” 
g^u- e = 1 is used to count powers of m, and labels such as the “n” in 
refer to perturbative order e”; I freely write these labels as either super- or 
subscripts. I work in geometric units with G = c = 1. 

1.1 A point particle picture: the MiSaTaQuWa equation and 
its interpretations 

1.1.1 Self-interaction with the tail of the perturbation 

As mentioned above, a point particle stress-energy tensor is not a well- 
defined source in the nonlinear Einstein equations. However, it is a fine 
source in the linearized theory. Eor simplicity, assume our object of mass 
m is isolated, such that in some large region, we can take it to be the 
sole source of stress-energy in the system. Now consider a metric = 
+ + where the background metric g^i, is a vacuum solution to 

the nonlinear Einstein equations, and is the leading-order perturbation 
due to our small object. Linearizing the Einstein equations in we obtain 

6Gf,,[h^]=87rT^^,, ( 1 ) 

where 6G^u[h] = I a-o linearized Einstein tensor, and is 

the leading-order approximation to the object’s stress-energy tensor. Now 
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suppose that the the object’s stress-energy can be approximated by 

T!^'^{x]z)= f mu^u'^6{x, z{T))dT, (2) 

J7 

which is the stress-energy of a point mass moving on a worldline 7 with coor¬ 
dinates z^ in the background spacetime. Here is the particle’s four- 

velocity, r is its proper time (as measured in g^v)i and d(x,z) = ^ 
is a covariant delta distribution, with g being the determinant of 5 ^ 1 /. 

Unlike the nonlinear Einstein equations with a point particle source, 
Eq. Q has a perfectly well-defined solution. If we introduce the trace- 
reversed perturbation — \g^vg'^^h\p and impose the Lorenz gauge 

condition = 0{e), then the linearized Einstein equation takes the 

form of a wave equation. 


(3) 

where 

E^,[h^] = nhl, + 2R^^/hip, (4) 

□ = g^'^V^Vu, and R^avp is the Riemann tensor of the background. As¬ 
suming the existence of a global retarded Green’s function G^j^y^iyi{x,x''^ 
for this wave equation, we can write the retarded solution as 

^ '' dV = A j dr, (5) 

where a primed index refers to the tangent space at the source point x'. 

Now, if the background were flat, waves would propagate precisely on 
the light cone, and the retarded Green’s function would be supported only 

on points connected by a null curve: ^^^ where 

(t, x“) is a Gartesian coordinate system, t — |x“ — x°‘'\ is the retarded time, 
and |x“ — = ^/dahix^^ — — x^') is the spatial distance between 

the points x and x'. However, in a curved spacetime, Huygen’s principle 
no longer holds. Waves scatter off the spacetime curvature, causing them 
to propagate from a source point x' both on the future null cone of x' and 
within that cone. Gorrespondingly, the Green’s function at a point x has 
support both on the past null cone of x and within it. This implies that 

^My conventions for the Green’s function are those of Ref. (3; Ref. also contains 
a pedagogical introduction to bitensors, objects which live in the tangent spaces of two 
different points x and x'. 
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if we look at the field at a point x near the worldline, we can split it 
into two pieces: a direct piece, corresponding to the portion of the field that 
propagated to x from a point z^(rret) along a null curve; and a so-called tail 
piece, corresponding to the portion of the field that propagated from 

all earlier points z^{t < Tret) on the worldline. The direct piece diverges 
like a Coulomb field, behaving as 1/r (where r is a geodesic distance from 
the particle). The tail piece, on the other hand, is hnite. 

A detailed analysis (such as the one presented in the bulk of this paper) 
reveals that at leading order, the mass m is constant, and the force on the 
particle vanishes; in other words, it behaves as a test particle in The 
same analysis applied at subleading order reveals that the direct piece of 
the field exerts no force on the particle, but the tail piece does exert a force, 
and the equation of motion is found to be 

^ (2'>S1 - Kt) + 0('"). (6) 

where ^ projects orthogonally to the worldline, 

and the tail term is given by 

hfupizir)) =Amf V u^'dr'; (7) 

J —OO 


the integral covers all of the worldline earlier than the point z^(t) at which 
the force is evaluated. The bar atop again denotes a trace reversal. 

Equation Q is termed the MiSaTaQuWa equation, after Mino, Sasaki, 
and Tanaka [23], who first derived it, and Quinn and Wald |26j . who re¬ 
derived it very shortly thereafter using a very different, independent method. 
The intuitive picture to glean from the MiSaTaQuWa result is that the direct 
piece of the field is analogous to a Coulomb field, moving with the particle 
and exerting no force on it, in the same way the self-held exerts no force on 
a body in Newtonian gravity. Very loosely speaking, from the perspective 
of the particle, the tail, consisting as it does of backscattered radiation, is 
indistinguishable from any other incoming radiation. In other loose words, 
it is effectively an external held, and like an external held, it exerts a force. 

Much of this paper is devoted to showing how the MiSaTaQuWa result 
can be robustly justihed, and higher-order corrections to it can be found, 
within a systematic expansion of the Einstein equations. As a byproduct of 
that analysis, in Sec. 4.1 I will show that the setup of the linearized system 
with a point particle source in this section is rigorously justihed—even for a 
non-material object such as a black hole, and even though at nonlinear orders 
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the field equations cannot be written in terms of such a source. However, 
before moving on to those matters, we may say significantly more on the 
basis of the point-particle result. 


1.1.2 The Detweiler-Whiting description: a generalized equiva¬ 
lence principle 


In the original MiSaTaQuWa papers [23l |26] , the authors noted that Eq. Q 
appears to have the form of the geodesic equation in a metric 
when that geodesic equation is expanded to linear order in However, 

is not in any way a nice field. It does not satisfy any particularly 
meaningful field equation, nor is it even differentiable at the particle [9] 
[despite superficial appearances in Eq. Q]. 

Detweiler and Whiting provided a more compelling form of the MiSa¬ 
TaQuWa result [MIES]. Rather than splitting the retarded field into a direct 
piece and a tail, they split it as 


h 


1 

fll/ 


-U h 


R1 


( 8 ) 


The singular field can be interpreted as the bound self-field of the parti¬ 
cle. Like the direct piece of the field, it exhibits a I/r, Coulomb divergence 
at the particle; but unlike the direct piece, it satisfies the inhomogeneous lin¬ 
earized Einstein equation = — IhvrT^j,, bolstering its interpretation 

as a self-field. Similarly, the regular field improves on the interpretation 
of the tail: it includes all the backscattered radiation in the tail, but it is a 
smooth solution to the homogeneous wave equation = 0. Hence, 

more than we could of the tail, we can think of as an effectively exter¬ 
nal field, propagating independently of the particle. Erom it we can define 
what I will variously call an effective metric or effectively external metric 

Fittingly, given this interpretation of h^y, Detweiler and Whiting showed 
that the MiSaTaQuWa equation (§ can be equivalently written as 

^ = - IfP"- + O(t^), (9) 


or (following Appendix [A]) explicitly as the geodesic equation in the metric 
9uv-> 

fS2 

,2\ 


df"^ 


= 0{e^), 


( 10 ) 


refer the reader to Appendix]^ for the expansion of the geodesic equation in powers 
of a metric perturbation. 
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where ^ ^ is a covariant derivative compatible with ^ 

is the four-velocity normalized in and f is the proper time along 
as measured in g^y. Equation Q is equivalent to Eq. Q because on the 
worldline, differs from h^^y only by (i) background Riemann terms that 
cancel in Eq. ([^ and (ii) terms proportional to the worldline’s acceleration, 
which can be treated as effectively higher order because the acceleration is 
already 

Allow me to dwell longer on the interpretation of the regular field. Be¬ 
cause g^y is a smooth vacuum solution, at the particle’s position an observer 
cannot distinguish it from g^y. Although a portion of g^j,y comes from the 
retarded field sourced by the particle, to the observer on the worldline, it 
appears just as would any metric sourced by a distant object. However, 
this interpretation of the effective metric as an effectively external metric is 
delicate. To effect the desired split into and both fields must be 
made acausal when evaluated off the worldline 13 [25]. More precisely, in 
addition to depending on the particle’s causal past, h^l{x) depends on the 
particle at spatially related points xb So in that sense its interpretability 
as a physical external field is limited. Yet when evaluated on the worldline, 
h^l and its derivatives are causal, and that is the sense in which g^y appears 
as a physical metric on the worldline. 

I impress upon the reader the significance of these properties of 
they are what makes Eq. (10) a meaningful result. Although it may not 
be an obvious fact at first glance, any equation of motion can be written 
as the geodesic equation in some smooth piece of the metric. This is most 
easily seen by writing the equation of motion in a frame that comoves with 
the particle. In locally Cartesian coordinates (t, x*) adapted to that frame, 
such as Fermi-Walker coordinates |9], the particle’s equation of motion reads 

is the particle’s covariant acceleration and is 


Qi = Fi, where = 


— 


~dF~ 


the force (per unit mass) acting on it. Now suppose we were to define some 
smooth field h^^y and a corresponding singular field = hj^y — K^i^y. In the 
comoving coordinates, the linearized geodesic equation in the regular metric 
+ h’^fjy, in the form analogous to Eq. ([^, reads a* = —h^li + 
matter what force F^ acts on the particle, the equation of motion o, = T) 
could be written as the geodesic equation in gfj,y + h'^^y simply by choosing 
= 0 and = 2Fi, for example. Besides those two conditions, 

the regular field could be entirely freely specified, and regardless of the 

specification we made, we would have defined a split h^y = h^^y + h 


rl 

''flU 


m 


which h^^y is singular and exerts no force, and gfj,y + h'^^jy is a regular metric 
in which the motion is geodesic. 
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Given this fact, it is of no special significance that the MiSaTaQuWa 
equation is equivalent to geodesic motion in some effective metric. But it is 
signihcant that the MiSaTaQuWa equation is equivalent to geodesic motion 
in the particular effective metric = g^i, + identified by Detweiler 
and Whiting, because of the particular, ‘physical’ properties of that metric: 
g^u is a smooth vacuum solution that is causal on the particle’s worldline. 
Because of those properties, we may think of the MiSaTaQuWa equation as 
a generalized equivalence principle: any object, if it is sufficiently compact 
and slowly spinning, regardless of its internal composition, falls freely in a 
gravitational field g^i, that can be thought of (loosely if not precisely) as the 
physical ‘external’ gravitational field at its ‘position’. I stress that this is a 
derived result, not an assumption, as will be shown in Sec.|^ By that stage 
in the paper, the principle’s reference to extended, “sufficiently compact and 
slowly spinning” objects will have become clear. 

Another aspect that will have become clearer is the non-uniqueness of the 
effective metric and the limitations of interpreting it in a strongly physical 
way. Nonetheless, the split of the metric into a self-field and an effective 
metric will be a recurring theme. In many ways, the generalized equivalence 
principle just described is both the central tool and the core result of self¬ 
force theory. It is conceptually more compelling than the expression for the 
force in terms of a tailit is less tied to any particular choice of gauge [19], 
and it is often easier to use as a starting point from which to derive formal 
results. Perhaps most importantly, it is easily carried to nonlinear orders: as 
I describe in later sections, at least through second order in e, the generalized 
equivalence principle stated above holds true naiiB]. 

1.2 Extended bodies and the trouble with point particles 

Although the point-particle picture seemingly works well within linearized 
theory, it is not obvious a priori how it fits within a systematic asymptotic 
expansion that goes to higher perturbative orders. For as I have reiterated 
above, one cannot use a point particle in the full nonlinear theory. Let me 
now expound on that point. 

® Another compelling physical interpretation is provided by Quinn and Wald [26]. They 
showed that the MiSaTaQuWa equation follows from the assumption that the net force 
is equal to an average over a certain “bare” force over a sphere around the particle. 
(This assumption was later proved to be true in a large class of gauges |20[|21[[55| .i In the 
language of Detweiler and Whiting fields, the force lines of the singular field are symmetric 
around the particle and vanish upon averaging, while the force lines of the regular field 
are asymmetric and add up to a net force on the particle. 
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Suppose we attempt to model the extended object as a point particle in 
the exact spacetime, with a stress-energy tensor 




dz^ — z°‘) , 

m— -;- ,_ — -dl, 


dt dt 




( 11 ) 


where t is proper time on z^ as measured in g^jy. If we expand the metric 
as g^v = g^u + Z)n>o ll^® linearized Einstein equation is exactly as 

described in the previous section. But at second order, severe problems arise. 
The second-order term in the Einstein equation reads 

h^] + ( 12 ) 


and is the second-order term in 


where 5‘^G^''[h,h] = 1 *^ 

. There are two problems with Eq. (12). Most obviously, contains 
terms like 




d^(x“ - z^) 


dr. 


As described in the preceding section, diverges as 1/r near the particle; 
hence, the stress-energy diverges in the distributionally ill-defined manner 
l(5(r). Even if we could somehow mollify the ill behavior of this piece of 
the source, the other source term in Eq. ( |12[ ) would still present a prob¬ 
lem. S‘^G^'^[h^, h^] behaves schematically as {dh^)‘^ + h^d^h^. Therefore, it 
diverges as 1/r^. Such a divergence is non-integrable, meaning it is well de¬ 
fined as a distribution only if it can be expressed as a linear operator acting 
on an integrable function. In the present case, there does not appear to be 
a unique way of so expressing it. 

These arguments make clear that a point particle poses increasingly 
worsening difficulties at nonlinear orders in perturbation theory. And it 
is well known that in any well-behaved space of functions there exists no 
solution to the original, fully nonlinear equation G^^[g] = SttT^^ with a 
point particle source mm- 

Despite these obstacles, one might suppose that the point-particle model 
could be adopted and the divergences resolved using post-hoc regularization 
methods. Evidence for this reasoning is given by the successful use of dimen¬ 
sional regularization in post-Newtonian theory [3]. In the fully relativistic 
context, the most promising route to such regularization appears to be ef¬ 
fective field theory jSH [35] . 

However, at a fundamental level, there should be no need for such reg¬ 
ularization in general relativity. Outside of curvature singularities inside 
black holes, everything in the problem should be perfectly hnite. For that 
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reason, in this paper I will be interested only in formalisms that deal with 
finite, well-defined quantities throughout. 

So let us do away with the fiction of a point particle and think of an 
extended object. Perhaps the most obvious route, at least at hrst glance, 
to determining this object’s motion is to go to the opposite extreme from 
the linearized point particle model, by looking instead at a generic extended 
object in fully nonlinear general relativity; if one is interested in the case of 
a small object, one can always examine an appropriate limit of one’s generic 
results. This line of attack has been most famously pursued by Dixon [36] 
(inspired by early work by Mathisson |37j) and Harte [38]. Specializing to 
material bodies, Dixon showed that all information in the stress-energy 
can be encoded in a set of multipole moments, and all information in the 
conservation law = 0 can be encoded in laws of motion for the 

body. These laws take the form of evolution equations for some suitable 
representative worldline in the object’s interior and for the object’s spin 
about that worldline. A “good” choice of representative worldline may be 
made by, for example, defining the object’s mass dipole moment relative 
to any given worldline and then choosing the worldline for which the mass 
dipole moment vanishes, establishing the worldline as a center of mass |39j . 
However, in order to transform these general results into practical equations 
of motion, Dixon’s method requires an assumption that the metric and its 
derivatives do not vary much in the body’s interior. Given that assumption, 
the force and torque appearing on the right-hand side of the equations of 
motion can be written as a simple expansion composed of couplings of the 
metric’s curvature to the object’s higher multipole moments; the higher 
moments, beginning with the quadrupole, may be freely specified, and their 
specification is entirely equivalent to a specihcation of the object’s stress- 
energy tensor. Equations of motion of this form can be viewed in Eqs. (13 )- 
(14) below. 

Unfortunately, for a reasonably compact, strongly gravitating body, the 
physical metric does not vary slowly in its interior, due to the body’s own 
contribution to the metric. To make progress, one must treat the object 
as a test body, an extended object that is non-gravitating but is equipped 
with a multipole structure—or, as discussed in Sec. 13 of Ref. [36|, one 
must hnd an efficacious means of separating the object’s self-held from the 
“external” held, analogous to the trivial split in Newtonian theory and to 
the Detweiler-Whiting split in the linearized point particle model. A well- 
chosen “external” held will vary little in the body’s interior, a well-chosen 
“self-held” will have minimal inhuence on the motion, and one can hope to 
arrive at equations of motion expressed in terms of couplings to the curvature 
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of the external field alone. 

In the fully nonlinear theory, finding such a split would seem to be highly 
nontrivial. Nevertheless, in a series of papers [281 EH iO] culminating in 
Ref. [38] (see also Harte’s contribution to these proceedings [30|), Harte has 
succeeded in finding a suitable split by directly generalizing the Detweiler- 
Whiting decomposition. He has shown that the “self-field” he defines mod- 
ihes the equation of motion only by shifting the values of the object’s mul¬ 
tipole moments, and the object behaves as a test body moving in the effec¬ 
tively external metric he defines. This extends the Detweiler-Whiting result 
from the linearized model to the fully nonlinear problem. However, there 
is one caveat to this generalization: beyond linear order, Harte’s effective 
metric loses one of the compelling properties of the Detweiler-Whiting field: 
it is not a solution to the vacuum Einstein equation. Despite this feature, 
Harte’s work is a tour de force in the problem of motion. 

The approach I take in this paper is complementary to Harte’s. Rather 
than beginning with the fully nonlinear problem, I will proceed directly to 
perturbation theory. There are several advantages to this. The perturbative 
approach naturally applies to black holes, while Harte’s formalism, because 
it is based on integrals over the object’s interior, is restricted to material 
bodies. The perturbative approach also naturally leads to a split into self- 
field and effective metric in which the effective metric satisfies the vacuum 
Einstein equation at all orders and is causal on the worldline. Most impor¬ 
tantly, the perturbative approach provides a practical means of solving the 
Einstein equations. In the fully nonlinear approach, one arrives at equations 
of motion given the metric, but not a practical way to hnd that metric. 

In the next section, I will begin to discuss the perturbative formalism. 
First, I note that with the work of Dixon and Harte, a new theme has been 
introduced: an object’s bulk motion can be expressed in terms of a set of 
multipole moments, and the i >2 moments are freely specifiable. Like the 
decomposition of the metric into a self-held and effectively external held, 
this second theme will appear prominently in the remainder of this paper. 

1.3 When perturbation theory fails near a submanifold: the 
method of matched asymptotic expansions 

As soon as we seek a perturbative description of the problem, we run into a 
new challenge. Say we assume an expansion of the exact spacetime (g^j,, Ale) 
about a background spacetime AIq); as in + O(e^), 

where all e-dependent terms are created by the small object (or by nonlinear 
interactions of its held with itself). This expansion assumes the object has 
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only a small effect on the metric. Clearly, if the object is compact, this 
cannot be true everywhere: sufficiently near the object, where r ~ e, the 
object’s own gravitational held will contain a Coulomb term ~ mjr ~ e®— 
just as large as the background and because it varies on the spatial 
scale e rather than e*^, having much stronger curvature than 

This fact motivates the use of matched asymptotic expansions. At dis¬ 
tances r ~ e® from the object (for example, r ~ M in an EMRI), we expand 
the exact spacetime around {g^u,M.Q), as above; I will call this the outer 
expansion. At distances r ~ e from the object (or r ~ m in dimensionful 
units), we introduce a second expansion, = g°^^ + + 0{e'^), where 

{9°uu ■, Mob]) is the spacetime of the object were it isolated, and the per¬ 
turbations are due to the helds of external objects (and to nonlinear 
interactions); I call this the inner expansion. In a buffer region around the 
object, dehned by e <C r <C we assume a matching condition is satished: 
if the outer expansion is re-expanded in the limit r <C and the inner ex¬ 
pansion is re-expanded in the limit r S> e, the two expansions must agree 
term by term in powers of r and e, since they both began as expansions of 
the same metric. The relationship between the various regions and expan¬ 
sions is shown schematically in Fig. and it will be described precisely in 
Sec. [21 

Historically, matched asymptotic expansions have been a highly success¬ 
ful way of treating singular perturbation problems in which the behavior of 
the solution rapidly changes in a localized region. For general discussions of 
singular perturbation theory in applied mathematics, I refer readers to the 
textbook |41j . and for more rigorous treatments, to Refs. |42[ I43j . For gen¬ 
eral discussions in the context of general relativity, I refer them to [Hill!. 

In the context of spacetimes containing small objects, matched expan¬ 
sions have been the standard method of tackling the problem; Refs. natta 
123 [22H211 are but a small sample. When it comes to obtaining equa¬ 

tions of motion, the method hearkens back to an early insight of Einstein 
and others lailMI]: an object’s equations of motion can be determined 
from the Einstein equations in a region outside the object. Specihcally, 
it can be determined from the field equations in the buffer region defined 
above. Relative to the object’s scale r ~ e, the buffer region e <C r <C is 
at asymptotically large distances, allowing one to make an asympotic char¬ 
acterization of how well “centered” the buffer region is around the object. 
One such characterization is based on a definition of multipole moments in 
the buffer region. Again, because the buffer region is at asymptotic infinity 
in the object’s metric g'^^, we can define the object’s multipole moments by 
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examining the form of the metric there, rather than having to refer to the 
object’s stress-energy. We can then use the metric’s mass dipole moment in 
the buffer region as a measure of centeredness: if we install a timelike curve 
7 in the background spacetime and define the mass dipole moment in coordi¬ 
nates centered on that worldline, then 7 is a good representative worldline if 
the mass dipole moment vanishes. Centeredness conditions along these lines 
will be described in more detail in Secs. and for a schematic preview, 
see Figs. and 

Prior to its application to self-force analyses [m m Eol this 

program was pursued furthest by Thorne and Hartle |48j . Where Dixon 
stood relative to the later non-perturbative work of Harte, Thorne and Har¬ 
tle stand in the same position relative to perturbative self-force construc¬ 
tions. As did Dixon, they derived general laws of motion and precession 
for compact objects. They considered an object immersed in some external 
spacetime, say and found forces and torques made up of couplings 

between the external curvature and the object’s multipole moments; specifi¬ 
cally, in Cartesian coordinates (t, x*) that are at rest relative to a geodesic 
of the external spacetime (and in which the mass dipole moment vanishes), 
they found 


— = + 0{e^ 


(13) 

(14) 


where and are the object’s linear and angular momentum relative 
to z^, and are its “mass and current” quadrupole moments (see 


the ends of Secs. 3.2 and 3.4), t is proper time (as measured in g^i,) on z^, 
and €ijk is the flat-space, Cartesian Levi-Civita tensor. The quantities (Bab 
and BBab are the electric-type and magnetic-type quadrupole tidal moments 
of the external universe, which describe the tidal environment the object 
is placed in; they are related to the Riemann tensor of q^u according to 
^ab = and BSab = -l€^'^{a^b)ipq- 

Physically, Eqs. (13)-(14) say that the object moves as a test body in 


the external metric^ the equations have the same structure as the test-body 
equations of motion mentioned above in the context of Dixon’s work j36| . 
Up to the coupling of the object’s moments to the tidal moments of the 
external universe, the motion is geodesic, and the spin parallel-propagated, 
in the external metric. But just as in Dixon’s work, these results become 
useful only once one knows how to split the full metric into a self-field and 
an effectively external metric. What Thorne and Hartle call the “external 
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metric” is not but rather it is essentially defined to be whatever creates 
tidal fields across the object. As Thorne and Hartle themselves remarked, 
this “external” field includes a contribution from the object’s own field. 
(Likewise, the “object’s” multipole moments can be altered by the external 
field.) 

In this sense, the self-force game is played by finding a useful split in 
which Thorne and Hartle’s general laws are valid. Doing so requires devel¬ 
oping a systematic theory of matched expansions for spacetimes containing 
small objects, as described in the body of this paper. Using those matched 
expansions, we will find that at linear order, we can circle back to the point 
particle picture: as first shown by D’Eath |45l I46j . the linearized perturba¬ 
tion in the outer expansion is identical to the linearized field sourced by 
a point particle (see also Refs. [HI [HI [20] for more refined derivations). 
And at that linearized level, we will find that the effectively external metric 

of Thorne and Hartle can be taken to be the Detweiler-Whiting effective 
metric = g^y + eh^l- 

1.4 Gauge, motion, and long-term dynamics 

Despite all the above preparation, we are still not ready to broach the prob¬ 
lem of motion in perturbation theory. Two additional facts must hrst be 
understood: in perturbation theory, motion is intimately related to gauge 
freedom [la EHEH; and in problems of astrophysical interest, the most 
important dynamical effects occur on the very long time scale 1 /e 153|. 

To the hrst point. At leading order, the object’s motion is geodesic in 
the background metric g^j^y] all deviation from that motion is driven by an 
order-e force. Suppose the self-accelerated worldline 7 is a smooth function 
of e. Then we can write its coordinates as an expansion 

z^^{s,e) = z^{s) + ez‘^{s) + 0{e^), (15) 

where s is a parameter on the worldline, and the zeroth-order term Zq is a 
geodesic of 5 ^ 1 /. Now consider the effect of a gauge transformation. Under 
a transformation generated by a vector a curve z^ is shifted to a curve 
z'^ = z^ — + O(e^). Nothing prevents us from choosing = z^, which 

leaves us with z'^ = Zq entirely eliminating the hrst-order deviation 

from Zq. This same idea can be carried to arbitrary order, meaning we can 
precisely set z'^ = Zq- In other words, the effect of the self-force appears to 
be pure gauge. 

In one sense, this result is true. If we look at any hnite region of space- 
time and consider the limit e —)• 0 in that region, the deviation from a 
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background geodesic is, indeed, pure gauge. This does not mean it is irrel¬ 
evant: in any given gauge, it must be accounted for to obtain the correct 
metric in that gauge. But it need not be accounted for in the linearized met¬ 
ric. We can always substitute the expansion (15) into to obtain 


ehl^ix; z) = ehl^{x; zq) -h e^Shl^{x; zq, zi) -h O(e^), 


(16) 


and we can then transfer the term 5hjj^^^{x; zq, zi) into the second-order 


perturbation, 
Eq. 


(An explicit expression for Shli^^{x] zq, zi) is given in 
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flU \ 


However, this analysis assumes we work in a fixed, finite domain—and 
as mentioned in the first paragraph of this section, we do not typically work 
in such a domain in problems of interest. Consider an EMRI. Gravitational 
waves carry away orbital energy from the EMRI at a rate E/E ^ e. It 
follows that the inspiral occurs on the time scale Er ~ E/E ~ 1/e, which 
is called the radiation-reaction time. So in practice, we are not looking at 
the limit e —)• 0 on a finite interval of time [0,T], where T is independent of 
e; instead, we are looking at the limit e —)• 0 on a time interval [0, T/e] that 
blows up. 

This consideration forces us to adjust our thinking about motion and 
gauge. Loosely speaking, the deviation from geodesic motion, ez^, is gov- 

On the radiation-reaction 
~ 1/e. In other words, it 


erned by an equation of the form ~ T/® 
time scale, it therefore behaves as ez^ 




blows up in the limit e —>■ 0. So on this domain, one cannot rightly write 
the worldline as a geodesic plus a self-forced correction, and one cannot use 
a small gauge transformation to shift the perturbed worldline onto a back¬ 
ground geodesic; the gauge transformation would have to blow up in the 
limit e —)• 0. 

Because we are solving partial differential equations and not ordinary 
ones, these arguments about time scales translate into arguments about 
spatial scales. For example, if we seek a solution in a Schwarzschild back¬ 
ground spacetime, fields (including inaccuracies in them) propagates out¬ 
ward toward null infinity along curves of constant t — r*, or inward toward 
the future horizon along curves of constant t -|- r*, where r* is the tortoise 
coordinate. If one’s accuracy is limited to a time span At, then it is also 
limited to a spatial region of similar size. 

To organize our thinking, let us denote by a spacetime region 

roughly of size ?(e) (both temporal and spatial). I call an asymptotic so¬ 
lution to the Einstein equations a “good” solution in if it is uniform 


^More precisely, it is governed by Eq. |35 
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in That is, the asymptotic expansion must 

eh^ 

satisfy lime_,.o = 0 and lime_s>o = 0 uniformly (e.g., in a sup 

C llnu 


norm). 

For the EMRI problem, we are interested in obtaining a good solution 
a domain Vi/^. Suppose we use an asymptotic expansion of the form (16) 


m 


and incorporate 5hfy{x] zo,zi) into In a gauge such as the Lorenz gauge. 


V 

and so 6h^ 


zf grows as 

"Di/g, its contribution to behaves at best as e 




^y{x-, zo, zi) likewise grows as Hence, on 
° comparable to gij,y. 


Clearly, this is not a good approximation. Suppose we instead eliminated 
zlf using a gauge transformation generated by = z^. This removes the 
offending growth in h‘j^y, but it commits a worse offense: it alters ehj^y by an 
amount which behaves at best as et, or as on Vi/^. Hence, if we 

are in a gauge where the self-force is nonvanishing, behaves poorly; if 
we are in a gauge where the self-force is vanishing, even behaves poorly. 

Let us chase the consequences of this. To obtain a good approximation 
in Pi/g, we need to work in a class of gauges compatible with uniformity in 
"Di/g. This means, in particular, that if we obtain a good approximation in a 
particular gauge—call it a good gauge —we must confine ourselves to a class 
of gauges related to the good gauge by uniformly small gauge transforma¬ 
tions. In turn, this means that the effects of the self-force are not pure gauge 
on Pi/g. Due to dissipation, z^ will deviate from any given geodesic Zq by 
a very large amount in Di/^, but by using an allowed gauge transformation 
we may shift it only by a very small amount, of order e, on that domain. In 
other words, although the self-forced deviation from Zq is pure gauge on a 
domain like T)^o, it is no longer pure gauge in the domain 


1.5 Self-consistent, Gralla-Wald, and oscnlating-geodesics ap¬ 
proximations 

In the preceding sections, we encountered several core concepts pertaining to 
the motion of a small object: (i) the metric can be usefully split into a self- 
field and an effectively external metric in which the object behaves as a test 
body, (ii) the bulk motion of the object can be described in terms of forces 
and torques generated by freely specifiable multipole moments, (iii) by using 
matched expansions, all of this can be done in a vacuum region outside the 
object, where multipole moments can be defined from the metric and laws 
of motion can be derived from the vacuum Einstein equation, and (iv) the 
representation of the object’s bulk motion, the asymptotic expansion of the 
metric, and notions of gauge must all be tailored to suit the long timescales 
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on which self-force effects accumulate. 

Let me now combine this information into a cohesive framework, mostly 
following Ref. [13] (but see also Refs. HMSHSH). The overarching method 
I describe consists of solving the Einstein equation with an “outer expan¬ 
sion” in some vacuum region outside the object, using only minimal 
information from the “inner expansion” to determine the behavior of the 
solution very near the object. 

Since we work in the vacuum region , we seek an asymptotic solution 
to the vacuum Einstein equation 


R^lv[^ — 0 - 


(17) 


Eollowing the lessons of post-Newtonian theory, I write this equation in a 
“relaxed” form 1 define impose the Lorenz gauge 

conditiorQ 

L,[h] = V^(V - = 0, (18) 

and write Eq. © with that condition imposed, making it read 


Ei,u[h] = Sfj,^[h], 


(19) 


where is the wave operator introduced in Eq. © (but here acting on 
the metric perturbation rather than on its trace reverse) and the “source” 
[g + h]— E^i, [h] is a nonlinear functional of . The background 
is chosen to be a smooth solution to R^u[g] = 0 (if matter exists outside 
the object, it is assumed to be sufficiently far away to lie outside This 

makes Eq. (19) a weakly nonlinear hyperbolic equation for the perturbation 
at this stage, that equation is still exact. 

I wish to solve Eq. (19) in subject to two types of boundary condi¬ 
tions: 


(BCl) Global boundary conditions. Examples of these are retarded bound¬ 
ary conditions or specihed Cauchy data. 

(BC2) The matching condition. In the buffer region, the solution must 
be compatible with an inner expansion. 

^Reference [HI describes how the entirety of this section can be performed in any gauge 
in which the linearized Einstein tensor is hyperbolic. Section]^ below offers a more general 
discussion of gauge. 

®In the Lorenz gauge in a vacuum background, the linearized curvature tensors are 
related by Efi,^[h] = —25Rfi,u[h] = —25Gfj,„[h] = Efj,u[h]. 
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Together these conditions ensure we are describing the correct physical sit¬ 
uation. But they do not yet uniquely determine the solution. Equation (19) 
is called “relaxed” because, unlike Eq. ( |17[ ), it can be solved no matter how 
the object moves; this relaxation arises because Eq. (19) is not constrained 
by the Bianchi identity, unlike Eq. The object’s motion is determined 
only once the gauge condition is also imposed, thereby making the solution 
to the relaxed equation also a solution to the unrelaxed one. As discussed 
earlier, in the present problem, the motion of the object is defined by the 
mass dipole moment of the metric in the buffer region, and we will find 
that the evolution of that mass dipole moment is determined by the gauge 
condition. 

Beyond these broad ideas, the specihcs of the method, and the size of 
the region in which it applies, depend crucially on how one repre¬ 

sents the object’s perturbed motion, which determines how one formulates 
the asymptotic solution to the Einstein equation. Here I describe three 
representations and corresponding asymptotic solutions: what I call “self- 
consistent”, “Gralla-Wald”, and “osculating geodesics” approximations. 


1.5.1 Self-consistent approximation 

In the self-consistent approximation, to avoid the secularly growing errors 


described in Sec. 1.4, I seek to directly determine an accelerated worldline 
Ze{s) that represents the object’s bulk motion; I do not wish to expand that 
worldline in powers of e. To accommodate this, I write the perturbation 

I /i, 

as hfj,iy{x,e', Ze), where z^ = and the quantities after the semicolon 
denote a functional dependence. I expand this functional as 


e; z,, Ze) = ^ i,). 


( 20 ) 


n>0 


Despite the fact that Ze and Ze depend on e, they are not expanded; in other 
words, I hold them fixed while taking the limit e —>• 0. Later, I will suppress 
the functional dependence on Ze and simply write h'^^{x]Ze). 

The self-consistent representation of motion is given its name because 
Ze must be determined simultaneously with /i”j^(x;Ze). It was the repre¬ 
sentation adopted in the original derivations of the MiSaTaQuWa equa¬ 
tion [23l [26l ES], and it is the one I used in describing that equation in 
Sec. o It was hrst put on a sound and robust basis, as part of a sys¬ 
tematic expansion of the Einstein equation, in Ref. [13]. In this section, I 
outline that expansion. 
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I first refine the region in which I seek a solution to the relaxed Einstein 


equation (19). Install the timelike curve 7 ^ (with coordinates z^) in the 
background spacetime, and let D 7 e be a region of size ?(e). I define 

to be a region of proper radius rs centered on 7 ^, with e <C <C 
The region I seek a solution in is then \ B^^. The inner 

boundary of this region lies in the buffer region, and there the solution must 
satisfy the matching condition (BC2). 


1.5.2 Field equations 


Since the coefficients h'^^{x-,Ze) in the expansion (20) depend on e, it may 
seem they are not uniquely determined. However, here I define the functional 


hnu{x,e;z) to be the solution to the relaxed Einstein equation (19); the 


relaxed equation places no constraints on z^, and each function 2 ;^ : M —)■ 
yields a different solution. (In the present context this means there are 
no constraints on the motion of the region B^.) The coefficients /i”^(x;z) 
are then uniquely determined to be the solution to the nth-order term in 
an ordinary power-series expansion of the relaxed Einstein equation. That 
nth-order term has the form ..., Up to n = 3, it 

reads 

Ef,t,[h^]=0 ( 21 a) 

E^i,[h^]=26‘^Rf,y[h^,h^] x G ( 21 b) 

E^^,[h^]=25^R^,[h\h\h^]+4S^R^^,[h\h^] ( 21 c) 

where I have defined the “nth-order Ricci tensor” to be the nth variatiorH 

1 d" 


6 Rfj,u[h, ..., /i] — ^^^R^iu[g + 


( 22 ) 


In the concrete calculations in this paper, I will require only 6‘^R^^, which 
is given explicitly by 

^ Ra/sih, h] = 2^^ T 2^^P’ hva^fi) 

T 4 ^^ \ahiiu\l3 2 ^^ hal3;^u ■ (23) 


Solving Eqs. (21) subject to the boundary conditions (BC1)-(BC2) yields 


a functional-valued asymptotic solution to the relaxed Einstein equation (19) 


Now we must find a particular function z^ = Ze for which h^u{x, e; z) is also 


^Cross terms like h^] are as defined in Eq. (227|. 


22 










an asymptotic solution to the unrelaxed equation (17). To do this, we must 
ensure h^y{x^ e; z) also satishes the gauge condition (18) to the same order 


I accomplish this in a systematic way by writing the accelerated equation 
of motion of Ze as 

" =F^{T,e), (24) 






and then assuming that like the force (per unit mass) appearing on the 
right-hand side can be expanded as 


F^(T,e) = 


n>0 


^F>i{T-,Ze,Ze 


(25) 


I substitute this expansion, together with the one in Eq. (20), into the gauge 
condition ( [I^ and solve order by order in e while holding (zt, Ze) fixed. By 
holding {ze,Ze) fixed during this procedure, rather than expanding their 
e dependence, I preserve the particular accelerated worldline that satisfies 
some appropriate mass-centeredness condition, such as the vanishing of a 
suitably defined mass dipole moment in the buffer region centered on 7 ^. 
Solving the sequence of gauge conditions to higher and higher order yields a 
better and better approximation to the equation of motion of that particular 
worldline, without ever expanding the worldline itself. The first few of these 
gauge conditions are 


L^[/i\Fo] = 0 , (26a) 

6 L^[h\Fi] = -L^[h^Fo], (26b) 

6 L^[h\F 2 ] = -L^[h^Fo]- 6 L^[h^Ffi- 6 ^L^[h\Fi,Ffi, (26c) 

where I have defined T] = and 

1 d^ 

6"L,lh,F,.. . ,F| = (27) 


This sequence determines the forces Ffii thereby determining the equation 
of motion (24). 


^Equation (2221 illustrates more explicitly, using a point-particle field, how the gauge 
condition implies an eqnation of motion. 

^^It also constrains other quantities in particularly determining the evolution of 
the object’s mass and spin. 
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1.5.3 Solution method 


We solve Eqs. (21) and (26) by working outward from the buffer region. 


Solving them in the buffer region using a local expansion, subject to (BC2), 
yields two things: the local form of the metric outside the object, and an 
equation of motion for the object. 

The local form of the metric is described in Sec. [HI In line with the 
themes of the earlier sections, it allows a natural split into an effectively 
external metric g^i, = g^i, + and a self-field where g^i, 

satisfies all the “nice” properties of the Detweiler-Whiting regular field. 

Derivations of the equations of motion at first and second order are 
sketched in Secs.[3|and[^ If at leading order the object’s spin and quadrupole 
moment both vanish, then the equation of motion is HSHH] 


^ - A" 


po-^u 


)uPu^ + 0(e3), 


(28) 


where + C>(e^). Following the steps of Appendix 






A.l 


this equation of motion can also be written as = 0 (e'*), the geodesic 

equation in the effective (smooth, vacuum) metric 5 ^,^. In other words, 
at least through second order in e, the generalized equivalence principle 
described in Sec. 11.1.21 holds. 

After obtaining the local results in the buffer region, one might think 


to solve Eqs. (21) and (26) globally in by imposing agreement with 

the local results on the inner boundary dB-f^ and then moving dB-y^ using 
the equation of motion. However, in practice, a global solution is instead 
obtained by analytically extending the buffer-region results into By^ \ 7 ^, 
replacing the physical metric there with the fictitious, analytically extended 
metric, while insisting that outside By^, the metric is unaltered. This pro¬ 
cedure (described in detail in Sec. allows us to work with field equations 
on the whole of U By^. At order e, the procedure reveals 

that z^) in is precisely equal to the perturbation produced by 

a point mass moving on Ze, as promised in Sec. 1 . 1.1 More generally, at 


all orders, it leads to a practical puncture scheme [I71I511I56H59], in which 
the puncture a local approximation to moves on 7 ^, and the field 
equations in By^ are recast as equations for a residual field that locally 
approximates 

'^^Note that although a puncture scheme utilizes approximations to and it is 
designed to exactly obtain (and any finite number of its derivatives) on the worldline, 
meaning it does not introduce any approximation into the motion of 7 e. Nor does it 
introduce approximations into the physical field h^i, = 
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The setup of a puncture scheme is compactly summarized in Eqs. (98)- 


(100) below. Using this scheme, one can directly solve for the effective 
metric on the worldline and use it to evolve that worldline via the equation 
of motion, and at the same time one can obtain the physical metric outside 
Although we begin with a potentially complicated extended object, 
this scheme illuminates the fact that in self-force theory, we do not need 
to know anything about the particularities of that object: at the end of the 
day, all necessary physical information about it is absorbed into the puncture 
and the motion of that puncture. 


1.5.4 Accuracy estimates 

How accurate will this self-consistent approximation be on a domain 
Let us make the reasonable assumption that the largest secularly growing 
error in the approximation arises from truncating the expansion ( |25[ ) at some 
order e"', leading to an error in of order 6z^ The largest 

error in h^i, is then 

Error in ^ ~ (29) 


The order of accuracy depends on the size of the domain we work in. Suppose 
we work in corresponding to the radiation-reaction time. On that 

domain, the error from neglecting is 0(e”'’'^T/e^) = Therefore, 

if we include only in the equation of motion, solving (via a puncture 
scheme) the coupled system comprising Eqs. (21a) and ([^, then the result 
contains errors 0{e^ • 1/e^) = 0{e )—which is as large as our first-order 
perturbation. In other words, this approximation fails on P ^/e- I f addition 
we include F^, solving the coupled system comprising Eq. (21a), (21b), and 
(28), then the error is 0{e^ ■ 1/e^) = O(e^); hence, with this approximation, 
we can have faith in our field /i^^(x; z). 

In the smaller domain corresponding to the so-called dephasing 

time tdph ~ the approximations are more accurate. Including only 

Ff, by solving Eqs. (21a) and ([^, yields a first-order-accurate solution with 
errors 0{e^ ■ 1/e) = c3(e^). Including F^, by solving Eqs. ( 21a[ ), (21b), and 
(28), yields a second-order-accurate solution with errors 0(e^ • 1/e) = 0{e^). 


'^^Here I return to what will become my common practice of dropping the subscript 
e on for simplicity, though I refer to the self-consistently determined center-of-mass 
worldline, not the freely specihable worldline for which the relaxed Einstein equation can 
be solved. 
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1.5.5 Gralla-Wald approximation 


I next consider the Gralla-Wald representation of perturbative motion, named 
after the authors of Refs. |20l I22j . In this approximation, as in Eqs. (15)- 


(16), one considers the effect of the self-force to be a small perturbation of 


the worldline, and one expands the worldline as 


‘(s,e) = ^eX(^), 


(30) 


n>0 


where the terms measure the deviation of from the zeroth-order 


worldline Zq. Substituting this expansion into Eq. (20), one obtains a new 
expansion of the metric perturbation: 


hf,^{x,e;z,z) = ^ e" 
n>0 


Ku{x]Zo, . . . ,Zn-l)- 


(31) 


Section [^describes this expansion in some detail. Explicitly, at first and sec¬ 
ond order, zq) = h^x; zq) and zq , zi) = ZQ)+6h]^^{x] zq, zi), 

as described below Eq. (1^). Since the individual terms Zn(s) are indepen¬ 


dent of e in Eq. (30), Eq. (31) is an ordinary expansion in which the coef¬ 


ficients do not depend on e, unlike in Eq. (|20|). Rather than starting 


from the self-consistent representation and then expanding the worldline, 
one could instead simply start with this ordinary expansion, as was done 
by Gralla and Wald. In this paper, to explicate the relationship between 
the two, I will instead almost always derive Gralla-Wald results from self- 
consistent results. 

In the Gralla-Wald approximation, the role of the domain is 

played by which excludes a small region around 70 (the worldline 

with coordinates Zq). In that region Eqs. (21) become 




0 X G 

2S^Rf,,[h\h\ h^] + A6^R^,[h\ h^] 


X £ V, 


70,<;(e)- 


(32a) 

(32b) 

(32c) 


The gauge conditions (26) become simply 


= 0 . 


(33) 
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The equation of motion (28) becomes a sequence of equations for zH ■ 


dTo 

dri 

D^4f 

dTo 


= 0 , 


= Fliro] Zo) - R'^alB-fUoZ^ul, 

= F^{to; zq, Zi) - R^a/3-f (uqZ^pU^ + 2ufz^u^'^ 


(34) 

(35) 


(36) 


dzH; 


Dz^ 


where tq is proper time on 70 as measured in Uq = and 

the forces and are constructed from according to Eqs. (206)- 

(207). The Riemann terms in these equations of motion are geodesic- 
deviation terms; they correspond to the fact that even in the absence of 
a force, two neighbouring curves z^ and Zq will deviate from one another 
due to the background curvature. Appendix |A .2 describes how Eqs. (34)~ 
(36) are derived from Eq. (28). As explained more thoroughly there, the 
quantities z^ and vectors that live on 79 . z^p is defined by apply¬ 

ing the expansion (30) in a normal coordinate system centered on 70 ; it is 
related to Z 2 in any other coordinate system by the coordinate-dependent 
relation zi^p = zl^ -|- \Typ{zQ)ziZp 

Just as in the self-consistent case, we can work with field equations on 
= ^ 7 o,<;(e) U B-yg by replacing the physical metric in B^yg with the ana¬ 
lytical extension of the buffer-region metric. Re-expanding the results from 
the self-consistent case, we find that zq) in 'D^g ,-^^^-^ is identical to the 

perturbation sourced by a point particle moving on Zg; the expansion of 
hU x; z) around zq) is derived in Sec. B We also arrive at a substan¬ 

tially simplified practical puncture scheme: rather than having to solve for z^ 
and hpy together, as a coupled system, one can first specify a geodesic Zq and 
then calculate in sequence (i) the perturbations hpi^{x;zo) and h^^{x;zo), 
(ii) the deviation Zy driven by (iii) the perturbations h‘j^,^{x-, zq,zi) and 
h^^{x] zq^zi), (iv) the deviation zi^, and so on. At all orders, the puncture 
moves on Zq; the deviations Zn> 0 ’ through their appearance in merely 

alter the singularity structure of the puncture. This puncture scheme is 
compactly summarized in Eqs. (117)-(119) below. 

On what domain is the Gralla-Wald approximation valid? There is no 
obvious estimate for the rate of growth of the terms but as in Sec. 


1.4 


we can easily estimate the growth of ezy to be of order If we assume 

^See Ref. [60] for an explicit solution to Eq. (351 in a particular scenario. 
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that the dominant error in /i^jy arises from the term in h?uy, then we have 


Error in h^i, 6 


21,1 






(37) 


If this estimate is valid, we can make the approximation valid to any order 
on domains with sufficiently small <^(e) <C tdph = l/v^- On domains 

comparable to the dephasing time, the deviation vector ez^ be¬ 

comes of order 1 , the second-order metric perturbation becomes as large as 
the first, and the expansion of the worldline ceases to be sensible. If higher- 
order deviations grow large much more quickly than z^, as we might 
surmise from Eq. (36), then the domain of validity of the expansion may be 


substantially smaller than even for first-order accuracy; if at each 

higher order the deviations grow more rapidly than the last, it may be the 
case that the domain of validity cannot be extended beyond 7^70,1 • 


1.5.6 Osculating geodesics 


Finally, I consider an approximation intermediate between the self-consistent 
and Gralla-Wald expansions, one which makes use of both the expanded 
and unexpanded representations of the worldline. Starting from the self- 
consistent representation, the idea is at each instant r on 7 , to perform 
a Gralla-Wald expansion in a region 117 ^^^^ , where 7o(t) is a geodesic 

of that is instantaneously tangential to 7 at time r; 7 o(t) is called an 


osculating geodesic [eng By solving the field equations of the Gralla-Wald 


approximation, one may calculate the self-force at time r, and then use that 
force to evolve z^ to the next time step. By following this procedure at each 
time step, one eventually obtains 7 over the entire timespan of interest. 

{x;z) 




The terms in the self-consistent approximation hpy{x]z) = 
can then be found simply by solving Eqs. (21a)-(21c) (and higher-order 
analogues) with 7 already pre-determined. 

More concretely, at each instant r on the worldline, one perform a Gralla- 
Wald expansion z^{T\e) = X] where r is the specific instant of 

interest and t' is variable, and one substitutes this expansion into the right- 


'^^The scheme I describe here should not be confused with the general method of oscu¬ 
lating geodesics, which is simply a way of using instantaneously tangential geodesics to 
rewrite an equation of motion " 7 ^ 7 !- = in terms of more convenient variables; that 
general method, inherited from celestial mechanics, is exact and does not inherently in¬ 
volve an expansion of , although it is particularly well suited to the osculating-geodesic 
approximation discussed here m- 
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hand side of Eq. (25 )—but not into the left-hand side —to obain 




— eF^{ t; Zo(^t)) + ^ -E^(r; 2 ;o(t)) 5^i(r)) + C>(e ). 


(38) 


The forces on the right-hand side are constructed from fields h^l{x] ZQf^^'j) 
and 2 ;o(t-), (and higher-order fields for higher order forces) ac¬ 

cording to Eqs. ( 206[ )-( 207| ). Eor example, in terms of the tail of the per¬ 
turbation, F{‘(r; Zo(r)) is given by the order-e term on the right-hand side 
of Eq. ([^ but with the tail integral ([^ evaluated over 


"0{r) 


rather than 


over z^. These fields are found by solving the following sequence of equa¬ 


tions at each value of r: (34) for the osculating geodesic z, 
h^ 




(x,Zo(t-)) in 


0(r)’ 


p5[ ) for zff^^y 


(37, . 2 o(r) 5 '^l(r)) 


m 


V. 


"yo(T) C^) * 


and (32b) for 

This suffices to compute F{‘(r; 2 ;o(t)) and E^(r; 2 :o(r), 2 i(t-)), but 


in principle the calculations could proceed to higher order. All of these equa¬ 
tions are to be solved subject to the “osculation conditions” = z>^{t), 
o(^)(r) = u^(r), = 0 , and = 0 , which state that 7 o(^) is 


u 


tangential to 7 at time r. They must also be solved subject to boundary 
conditions that ensure the sum X)n>o ^o(t)) agrees with the full field 

hfj_i, in at nonlinear orders, finding those boundary conditions may 

be highly nontrivial. 

At linear order, the osculating-geodesic approximation has already been 
concretely implemented to find z^ using the force Ef (r; Zo(t-)) [fi2ti64] (build¬ 


ing on the framework in Ref. 


and to compute h}^^{x; z 


How¬ 


ever, to my knowledge, the brief sketch above is the first time it has been 
described at nonlinear orders (although possibly equivalent ideas have been 
presented by Mino [HS])- Considerably more work must be done to estab¬ 
lish that the scheme is viable beyond linear order. If it is, then one may 
naively estimate that it is valid to the same order on the same domain as the 


self-consistent expansion: including only {t; Zq(^t-)) in Eq. (38) yields an 
approximation valid up to C>(e^) errors in a domain and including 

both F^{t; Zo(t)) and E^(r; 2 :o(r)) H(r)) yields an approximation valid up to 
0{e‘^) errors in a domain 


1.6 Outline of this paper 

As discussed above, the bulk of this paper focuses on the formalism presented 
in Ref. [13] and further developed in the series of papers [HI [TBHIHl El] • 

I begin in Sec. with a more complete description of matched asymp¬ 
totic expansions, the concepts of which underly most of self-force theory. I 


29 















focus on the importance of the buffer region and formulating appropriate 
definitions of the small object’s representative worldline. 

Sections and present an algorithm for constructing an nth-order 
expansion in the self-consistent approximation. In Sec. I describe how to 
obtain the general solution in the buffer region. This algorithm determines 
the equation of motion as well as a natural split of the general solution into 
a self-field and an effectively external field. In Sec. I describe how to 
generate the global solution in using as input the self-field obtained 

in the buffer region. Along the way, I show how at linear order, the point- 
particle picture is recovered, and how at nonlinear orders, a certain point- 
particle “skeleton” can be defined to characterize the object’s multipole 
structure. 

In Sec.[^ I show how to recover a Gralla-Wald expansion or an osculating- 
geodesics expansion from the results of the self-consistent expansion. 

In Sec. I discuss the gauge freedom in each of the three types of 
expansions—self-consistent, Gralla-Wald, and osculating-geodesics—and how 
that freedom relates to the representation of perturbative motion in each. 

In Sec. I sketch how one can derive equations of motion by obtain¬ 
ing further information about the inner expansion and then utilizing the 
relationship between gauge and motion. This is an alternative to the algo¬ 
rithmic approach in the buffer region, and it is the method that has been 
used in practice to derive second-order equations of motion UniESlES!. 

I conclude in Sec. with a summary and a discussion of future directions. 

The appendices contain more general results: expansions of the geodesic 
equation in a perturbed spacetime, expansions of point particle functionals 
of an accelerated worldline around a geodesic, and some identities pertaining 
to gauge transformations of curvature tensors. 

Some portions of this paper first appeared in slightly different forms in 
Ref. [H]: specifically. Appendices [A| [B| and and parts of Sec.j^ 

2 Matched asymptotic expansions 

2.1 Outer, inner, and buffer expansions 

Traditionally in applied mathematics, the method of matched asymptotic 
expansions has been used to find both inner and outer expansions to some 
desired degree of accuracy and then combine them to obtain a uniformly 
accurate solution in the entire domain. But for my purposes here, I will 
not be interested in obtaining an accurate solution in the region r ~ e; 
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instead, I will be interested in the inner expansion only insofar as it con¬ 
strains the outer expansion. Finding an accurate inner expansion would 
require specifying the type of compact object we are examining, while here 
I am interested in generic results that apply for any compact object. As 
it turns out, finding those generic results requires only minimal knowledge 
of the inner expansion—in fact, just its general form in the buffer region 
suffices. 

To build toward that conclusion, let me describe the formalism more 
geometrically. I focus for the moment on the self-consistent case. In the 
outer expansion, we wish to approximate the exact spacetime (g^i,,Ade) in 
the domain outside a small region of size r <C around the object. 

I expand the exact spacetime around a background by adopting 

some identification (p : C AIq —^ -M-e between the spacetimes (in the 

region outside the inner region). In a given coordinate system : Aio ^ 
the identification map assigns points p G Aio and ip{p) G the same 
coordinate values x^{p). I next write as a functional of a worldline 
7 C Aio as g^iy(x,e;z) in those coordinates, where = x^( 7 (s)), and 

I then expand for small e while holding both x^ and fixed. I thence 
arrive at Eq. (20),^ Note that is defined in the background, not in 
the perturbed spacetime; the identification map is assumed to exist only 
in not in the inner region where there may lie a black hole rather 

than an identifiable worldline, for example. Later, I will replace the physical 
metric in the interior of the object with the effective metric g^u, and via the 
identification map, the worldline will have identical coordinate values in the 
effective spacetime {g^i,,Ai) as in the background. 

Now, the inner expansion is constructed by choosing some coordinates 
(t, X*) centered on 7 , and then rescaling spatial distances according to x* = 
X*/ e. The inner expansion is performed by expanding for e —?• 0 while holding 
the scaled coordinates {t, x*) fixed, as i-CH 

g„„(i,S“,e) (39) 

n>l 


Tensors in this expansion live on a manifold Adobj; where Adobj) is 

truth, it is unlikely that any of the expansions I consider, whether self-consistent, 
Gralla-Wald, or osculating-geodesics, is convergent. More likely, they are asymptotic 
approximations only. So when performing the self-consistent expansion, I actually assume 
that |gj,,.(a;,e)-g))),(a;,e; 2 ;)| = o(e^), where e; z) = g^L,,{x) + J2n=i z). The 

notation o(fc(e)) means “goes to zero faster than fc(e)”. 

^^Here indices refer to the unsealed coordinates If components are written in 

the scaled coordinates, overall factors e and appear in front of ta and ab components, 
respectively. These overall factors have no practical impact. 
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identified as the object’s spacetime were it isolated. The perturbations 
z) describe the effect of interaction with the external spacetime. 
What is the meaning of this expansion? The scaled coordinates serve to 
keep distances fixed relative to the object’s mass in the limit e <C 1, effec¬ 
tively zooming in on the object by sending all distances much larger than 
the mass off toward infinity. The use of a single scaling factor makes the ap¬ 
proximation most appropriate for compact objects, whose linear dimension 
is comparable to their mass. Scaling only distances, not t, is equivalent to 
assuming the object possesses no fast internal dynamics; that is, there is no 
evolution on the short timescale of the object’s mass and size. 

Note that the treatment of the “region around the object” depends 
strongly on whether one considers a self-consistent or Gralla-Wald expan¬ 
sion. In the self-consistent case, the outer expansion takes the limit as 
the object shrinks toward zero size around the self-consistently determined, 
accelerated worldline 7 C Wfoj and the inner expansion blows up a region 
around that accelerated worldline; in the Gralla-Wald case, the outer expan¬ 
sion takes the limit as the object shrinks to zero size around the zeroth-order, 
background geodesic 70 , and the inner expansion blows up a region around 
that background geodesic. 

In either case, the relationship between the two expansions is illustrated 
in Fig. I now use this relationship to feed information from the inner 
expansion out to the outer expansion. This exchange of information is done 
in the buffer region around the object. From the perspective of the inner 
expansion, the buffer region lies at asymptotic infinity in Alobj- In that 
region, the inner expansion can be expressed in unsealed coordinates as 

gi,u{t,x°‘/e,e) (40) 

n>l 

and then re-expanded for small e (or equivalently, expanded for r ^ e; i.e., 
for distances that are large on the scale of the inner expansion). Gonversely, 
from the perspective of the outer expansion, the buffer region lies in a tiny 
region around the worldline. Hence, in that region the outer expansion can 
be expanded for r <C (i-e., for distances that are small on the scale of 
the outer expansion). Since the inner and outer expansions are assumed 
to approximate the same metric, it is assumed that the results of these re¬ 
expansions in the buffer region must match order by order in both r and 


'^^This matching condition amounts to the assumption that nothing too “funny” happens 
in the buffer region. It can instead be replaced by more explicit assumptions on the 
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Now, say the nth-order outer perturbation is expanded as = 

YlprPh^f,(t,n^) in the buffer region, where n* = x^jr.^^ One must allow 
negative powers of r, since at least part of the field will fall off with distance 
from the body. But there is a bound on the most negative power at a given 
order in e: Since the inner expansion is assumed to be well behaved, it 
must include no negative powers of e. And since any term in must 

coprrespond to a term in the inner expansion, if is written as a function 

of the scaled distance f = r/e, it must likewise have no negative powers of 
e. From this it follows that 




—n 


(t,n*) + 0(eV 




(41) 


Any higher power of 1/r would induce illegal powers of e; for example, 
1 


K+T 




IS 


m 

n,—n 


We can go one step further with this general analysis. Since e”/r 
independent of e in the scaled coordinates, and is the only term 
the inner expansion (39) that does not depend on e, it must be that h^u 
is equal to a term in If we write g°^\t^x^) in terms of the unsealed 

coordinates and expand for r ;§> e, we find its form in the buffer region ia^ 


5°^j(t,xV6) = 

n>0 ' 


(42) 


From the matching condition, we then have 

un,—n „obj,n 
= gf,u . 




(43) 


Therefore, at each order in e, the most singular (as a function of r) piece 
of the metric perturbation in the buffer region is determined by the 
r ^ e asymptotic behavior of the object’s unperturbed metric. In addition, 
note that because quantities in the inner expansion vary slowly in time 
relative to their variation in space, they are quasistationary: after changing 
to scaled coordinates, a derivative with respect to t effectively increases 


behavior of the full metric g^i^, such as the conditions assumed in Ref. or various 
others discussed in Ref. |43 |. 

^®lnr terms also generically arise. For simplicity, I incorporate those terms into for 
the moment. Their presence does not spoil the well-orderedness of the expansion, since 
r^llnr)"^ ^ (Inr)"* for p > p'. Similarly, Ine terms can occur in solving the relaxed 

Einstein equation and I absorb them into the coefficients z). 

^°The fact that the inner background must be asymptotically flat, containing no positive 
powers of r, follows from the assumption that the outer expansion contains no negative 


powers of e, in the same manner as the cutoff on powers of 1/r in Eq. (411. 
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Figure 1; The exact spacetime (top left), external background spacetime 
(top right), and inner background spacetime (bottom), with maps between 
them. In the exact spacetime, we have some small compact object, here 
shown as a material body (in dark red and black), but possibly a black 
hole or more exotic object. Around the object is a buffer region, shown in 
orange. Points outside the object can be identified (via the identification 
map (^g) with points in the external background manifold A^o- In a chart 
(t, X*) : A^o — ^ centered on a worldline 7 G A^Oj the points p G AIo 

and ipe{p) £ A4e are assigned the same coordinate values (t(p), x*(p)). The 
map ° zooms in on the object by rescaling spatial coordinates; in the 
same chart (t, x*), the point q = is assigned coordinate values 

{i{q),x^{q)) = (t(p),x*(p)/e). 


the power of e relative to a derivative with respect to x*. Hence, on any 
short time, we can think of being stationary, and we can write its 
asymptotic form (42) in terms of a canonical set of multipole moments |66l 
ET]. So our final statement is that the inner expansion constrains the outer 
expansion to have the form (41), and the coefficients can be expressed 

in terms of multipole moments of the small object’s unperturbed spacetime 


Adobj)- This is all the information that will be required from the i 
expansion (until we get to Sec. [^. 


inner 
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2.2 Defining the worldline 


Let us now use the above formalism to define what we mean by the object’s 
worldline. First consider the self-consistent case. In the coordinates (t, x*) 
centered on 7, calculate the mass dipole moment M* of the spacetime 
A mass dipole moment indicates the position of the center of mass relative 
to the origin of the coordinates. If the coordinates are mass-centered, a 
Coulomb-like piece of the field behaves as m/r; if the origin lies slightly 
away from the center of mass, by an amount then the Coulomb-like field 
behaves as m/|x® — Expanding this around = 0 , we find mjr + 
-|-O(I^P). The quantity is the mass dipole moment ML Ergo, 
if this mass dipole moment vanishes in the coordinates centered on 7, then 
the object is appropriately centered “on” 7, and we identify 7 as a good 
representative worldlineBecause the definition only utilizes quantities in 
the buffer region, this definition makes sense even if the object is a black hole 
or contains topological oddities such as a wormhole: even if there exists no 
identification map between the background and the exact spacetime in the 
region inside the object, at least in the buffer region the coordinates (t, x*) 
can be used to chart the manifolds Alg and Aioh]-, and in those coordinates 
the metrics and have reference to the worldline 7 that is defined 
only in AIq- I refer the reader again to Fig. [^to illuminate this. 

Now consider the Gralla-Wald case. Here 70 is the worldline around 
which the inner expansion is performed. In a generic gauge (in particular, 
in the Lorenz gauge I work in), the small object’s center of mass deviates 
from this worldline, as assumed in the expansion ( 30 ). Gauges in which 


the object does not deviate from 70 (on short timescales) are discussed in 
Sec.[^ but in a generic gauge, clearly we do not have M* = 0 in coordinates 
centered on 79. Instead, we calculate M* in those coordinates and then 
define the first-order correction to the motion to be 


z\ = My-m. ( 44 ) 

Figure [^illustrates the difference between this setup and the self-consistent 
one. 

^^This notion of mass-centeredness based on the mass dipole moment of applies only 
to order-e deviations from z^. For higher-order deviations, mass-dipole-moment terms in 
the perturbations must also be considered, or some other copacetic centeredness 
condition must be imposed, as discussed in Secs. [Sj and Pfj 
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Figure 2: Perturbative motion of a small object as defined in the Gralla- 
Wald (left) and self-consistent (right) approximations. In the Gralla-Wald 
case, the buffer region is centered on a zeroth-order worldline 70 C Mq, 
which in a given coordinate system : AIq 1^^ has coordinates Zq(s). In 
a generic gauge, the self-force drives the object (shown as a sequence of red 
circles) away from 70 , and this deviation is represented by Zi{s), a vector 
defined on 70 . In the self-consistent case, the buffer region is centered on 
an accelerated worldline 7^ C Mo, which in the same coordinates has 
coordinates z^{s, e). The accelerated worldline faithfully tracks the object’s 
motion, such that the object is always at the “center” of the region enclosed 
by the buffer. On timescales much smaller than the dephasing time, we have 
that the two approximations are related by z^{s, e) = ZQ{s) + €z^{s) + 0{e‘^). 

3 Algorithm for an nth-order self-consistent ap¬ 
proximation: general solution in the buffer re¬ 
gion 

We are now positioned to actually obtain an outer expansion. In this section, 
I present an algorithm for finding the outer expansion in the buffer region. 
In Sec. 1^ I describe an algorithm for obtaining a global solution using as 
input the results from the buffer region. In both cases, I specialize to the 
self-consistent case. 

My method of finding the general solution in the buffer region is mod¬ 
eled on the post-Minkowskian methods of Blanchet and Damour [ 68 ] . Like 
them, I write the general solution in terms of a set of algorithmic multipole 
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moments, which are defined simply from the algorithm of solving the relaxed 
Einstein equation, but a subset of which can be concretely identified with the 
physical multipole moments of . All the moments are completely uncon¬ 
strained so long as we solve only the relaxed Einstein equations (21). Once 
we impose the gauge condition, as in Eqs. (26), the moments become con¬ 
strained, and in particular, evolution equations arise for the mass monopole, 
mass dipole, and spin dipole moments. The equation for the mass dipole 
moment will be used to identify the center-of-mass worldline 7 for which 
M* = 0 . 

The solution thus obtained will be a general solution in the sense that it is 
made to satisfy the matching condition (BC2) but not any particular global 
boundary conditions (BCl). (Refer back to the opening portion of Sec. 1.5 ) 
Because the explicit calculations, as well as their results, are exceedingly 
lengthy, I merely sketch the algorithm and the form of the results. I refer 
the reader to Refs. p la im Eu for more detailed expositions. 


3.1 Setup 

Eor concreteness, I work in Fermi-Walker coordinates {t, x°‘) centered on 7 , 
in which t is proper time on 7 , r = is the proper distance from 

7 along a spatial geodesic (3 that intersects 7 perpendicularly, n* is a unit 
radial vector that labels the direction along which /3 is sent out, and the 
spatial coordinate is x* = rn®. Ref. [9] contains a pedagogical introduction. 


The background metric in these coordinates is given by 

gtt = -(1 + ttix'- f - RoiQjx'^x^ + 0{r^), (45a) 

gta = -iRoiajX^’x^ + O(r^), (45b) 

9ab = Sab - iRaibjx'x^ O(r^), (45c) 

where the Riemann terms are evaluated on the worldline and contracted 
with members of a tetrad (u",e") on 7 that satisfies for any 


coordinates y^. For example, Roiaj{t) = Ra^isy{z{t))u'^ete't. An overdot 
will indicate a covariant derivative along the worldline, such as Roiaj = 
Ra^j.i 3 u;p\^u^&i e-aS-’jU'^■ Relating back to the discussion below Eqs. 
the Riemann tensor on 7 can be written in terms of two tidal moments: 

^ab — RoaObj ^ab — {a^b)Qpqi (^^) 

and Rabcd ^acSbd T ^bd^ac T ^ad^bc T ^bc^ad' 


(13)-(14), 
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Reference m displays the metric ( |45[ ) to higher order in r. At all orders 
in r, all terms are made from the acceleration, Riemann tensor, and deriva¬ 
tives of the Riemann tensor. Because the coordinates are tethered to an 
e-dependent worldline, the background picks up an e dependence, seen most 
obviously in the acceleration terms a*. However, this is purely a coordinate 
effect. In any global coordinates covering V^^c;{e)i 9^lu is e-independent by 
definition. In a Gralla-Wald expansion, the e-dependent coordinate trans¬ 
formation to Fermi-Walker coordinates would be expanded in powers of e, 
splitting it into an e-independent coordinate transformation (to Fermi coor¬ 
dinates centered on 70 ) plus a gauge transformation. Because I wish to at 
no point expand 2 ;^, I do not expand the transformation in this way. 

In these coordinates, I expand as in Eq. (41). I further assume 
that the coefficients in that expansion are smooth, and I make logarithms 
explicit; these logarithms appear generically in solutions to inhomogeneous 
hyperbolic equations We then have an expansion of the form 

^max 9max 

hly{t,x'^-z)= ^ ^ rP(lnr)5/i"^f(t;z)h-^, 

p>—n £=0 g=0 


(47) 


where L = ii ■■■ i£ is a multi-index, is a smooth function of t, and 

is an STF combination of unit vectors. The decomposition 
in terms of is equivalent to an expansion in spherical harmonics, and 
like a spherical harmonic, is an eigenfunction of the flat-space Laplacian, 
satisfying d^dth^ = — References | 68 l I70j contain excellent intro¬ 

ductions to this type of decomposition, along with many useful identities. 
Favoring n“ over angles and over spherical harmonics, is useful 

because is naturally written, as above, in terms of = rre“. 


3.2 Seed solutions 

By definition, is the solution to the nth-order relaxed Einstein equa¬ 
tion, which has the form ,... Before considering 

^^Intuitively, the logarithms are caused by the object perturbing the spacetime’s light 
cones. One can expect the solution to the exact Einstein equation to propagate on (and 
within) null cones of the exact spacetime, and given that the mass of the body induces a 
logarithmic correction to the retarded time, logarithmic corrections then naturally appear 
in This effect is well known from solutions to the Einstein equation in harmonic 

coordinates (see, e.g.. Refs. [681169] '). For generality, I allow logarithms at any value of 
n, but I assume that for each finite n, p, and £, the highest power of Inr is a finite 
number qniii.x{n,p, £). For simplicity, to make sure that term-by-term differentiation is 
valid without worrying about issues of convergence, I also assume for a given, finite n and 
p, £ has a maximum £ina,x{n,p). 
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the general solution to this equation, let us first consider solutions to the 
corresponding homogeneous equation Efj,^[h'^] = 0. We shall find that every 
solution of the form (47) is a member of a set z,I^) : 0 < £ < n} 

or of a set (x; z, fc”) : £ > 0 }, or it is a superposition of such members. 
For a given background metric, the seed solutions z, /”) and the free 

solutions are completely determined by the worldline and 

functions of time or s-long that worldline. The quantities 

identified as algorithmic multipole moments that describe the 
object’s multipole structure, while each of the quantities will identify 
an effectively external, freely propagating perturbation. 

We arrive at these conclusions simply by substituting the expansion (47) 
into = 0. First note that spatial derivatives reduce the order in r 

by one power, while temporal derivatives do not affect the order in r. This 
means we can write Efj,,y[h'^] = d'^dih'^^ + where ~ h^^/r. 

Solving for order by order in r is hence reduced to solving a sequence 
of Poisson equations d'^dih’^^, = — Substituting the expansion (47) 

into this, we find that the coefficient of rP“^(lnr)'^h^ in the equation is 


b(p +1) - e(t + + (9 + l)(2p + + (9 + 1)(9 + 

p-2,q,e 


< n,p,q+l,£ 


n,p,q+2,i 


P 1 ^max ^max 

E EE'^"^. 

p'=—n 


r'P (Inr)*^ n 


L' i^np'q't' 
'‘‘pvL' 


(48) 


where means “pick off the coefficient of r^“^(ln rYh^". The important 

thing to note is that the right-hand side depends only on coefficients with 
p' < p. So start with the term with the lowest power of r, such that 

the right-hand side of Eq. (48) vanishes. Further, start with q = Qraax, and 
suppose gmax > 0. By assumption, ^n^-n.,gmax+i,^ vanish, 

and so Eq. (48) becomes [n(n — 1) — £{i + = 0. This has a 

nontrivial solution only if £ = n — 1. Now look at the equation for the same I 
and p but for g = gmax-1- It reads [n(n-l)-£(£-t-l)]h”;;2"’'^“""~^’^+9max(l- 
~ term vanishes because I = n — Hence, 

= 0 if ( 7 max > 0 , and we conclude that is 

the only nontrivial solution to the 1 term in the homogeneous equaion 
Efj,u[h'^] = 0. I define the algorithmic moment 


rn _ ,n,—n,0,n—l 
^puL — '^pi/L 


(49) 


Now proceed to sequentially higher orders in r, and at every order, set to 
zero all free functions that arise. The result is a solution to E^uih'^] = 0 of 
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the form 




rn 

^fiuN 




N-l 






(50) 


in which every single term is linear in (and t-derivatives of it), 

specifically consisting of (and its derivatives) contracted with the 

other available Cartesian tensors n*, a*, (5®-^, pieces of the background 
Riemann tensor, and pieces of derivatives of the background Riemann tensor. 

Now set to zero and solve the l/r®®+^ and higher-order terms in 

Eq. (48), following exactly the same steps as above. Again the result is a 


homogeneous solution that is completely determined by a single function, 
in this case I'^yjy [_2 = ■ Now keep doing the same by beginning 

at each following order up to (and including) 1/r^, always setting every free 
function to zero except the first one encountered. This procedure completely 
populates the set {h^^y'^{x; z, I^) : 0 < £ < n} with solutions of the form 


Lseed/ . Tn\ _ > 

^LIU ) D. 




+ 0 ( 1 // 


(51) 


where every term in the solution is proportional to 

or its derivatives. The algorithmic multipole moments are symmetric 
in their first two indices and STF in their last £ indices. 

Continuing the same procedure to higher orders in r, we completely 

,free 


populate the other set, z, I”) : £ > 0}, with solutions 


hUr (^; + fc") = r%yL{t)n^ + 0 (/+/ 


Juri 




(52) 


where every term in the solution is proportional to kUit) = or 

its derivatives. The change in behavior from to arises because the 


factor \p{p + 1) — £{£ + 1)] in Eq. (48) vanishes for p = —— lifp<0 and 


for p = £ if p > 0. As with the quantities are symmetric in their 
hrst two indices and STF in their last £ indices. 

It is easy to see from the above steps that the set of seed solutions 
z, I”) and free solutions h^(x; z, k^) form a complete basis of solu¬ 
tions to E^y[hr‘] = 0 of the form ( |47| ). If extended down to r = 0, each seed 
solution diverges at r = 0, while each free solution is smooth there. 

For later purposes, it will be convenient to split each algorithmic multi¬ 
pole moment into a mass and current moment, 


jn _ T\/rn , an 


(53) 
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where the mass moment is the even-parity part of satisfy¬ 
ing the current moment is the odd-parity 

part, satisfying S'^ji = This splits each seed solution into two, 


Kf{x-z,M^) and Kf{x-z,S^). 


3.3 General solution 

Now turn to the inhomogeneous equation E^,y[h'^] = ..., Start 

at the n = 1 equation, (21a). The source vanishes, and the general solution 


is made up of a single seed solution, 
W^f{x-zjl) = 




+ 0(r°), 


plus the sum of all free solutions z, k^). That is, 

hU = 4') + E ^1)- 


(54) 


(55) 


£>0 


Next move to the n = 2 equation, (21b). The l/r^'*'^ term in the equa¬ 
tion looks essentially the same as Eq. (48), except that the right-hand side 
contains a term coming from [h }, One can straightforwardly solve 
for the functions to find that each one of them is either the starting 


point for a new seed or free solution or for £ > 0) 

or it is directly proportional to a pqi mode of the source. Since the source is 
constructed from quadratic combinations of it follows that the general 
solution is made up of the two seed solutions 


h^^f{x;z,li) 



+ 0{l/r), 



+ O(r 0 ), 


(56) 

(57) 


plus the sum of all free solutions /i|)'®®(x; z, kj), plus a particular inhomoge¬ 
neous solution made up entirely of terms that are quadratic combinations of 
members of the set ^jiuL • ^ ^ 0} (™ other words, a particular solution 
in which all seed solutions and free solutions are set to zero). That is. 


h 


2 

flU 


h^^fix; z, If) + h^^fix- z, II) + E z, kj) 

e>Q 

+ hli^{x-,z,I^,{k}}i>o). 


(58) 
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The promised logarithms of r first appear in {x\ z, Iq, {kj}e>o) [IS [IT]. 
Carrying this procedure to nth order, we find 

n—1 

h% = E -'?) + E *■•?) 

f=0 £>0 

+ {Ii^}m<n,l<n, {k^}m<n,e>o), (59) 

where all functions appearing in h'^i^ are made up of explicit products 
of members of the set : m < n, £ < n — 1} U {k'^^i^ : m < n, f > 0} (or 

their t derivatives) contracted with combinations of n*, a*, 5*-^, pieces of 
the background Riemann tensor, and pieces of derivatives of the background 
Riemann tensor. 


In any particular solution of the form (47) to the relaxed Einstein equa 


tion outside a small object, every term, to all orders in r and e, is a lin¬ 
ear or nonlinear combination of the functions 

Specifying all the functions /"^.^(t) and kl^^^{t) corre¬ 
sponds to making a specific choice of small object and of global boundary 
conditions. 


3.4 Imposing matching and gauge conditions 

To satisfy the unrelaxed Einstein equation, the general solution described 


above must also satisfy the gauge conditions (26), which we can write gener- 
ically as 

L,[h\F,] = P,[{h^,F,^r^},], (60) 

for some that is linear in the /I’s but not usually in the P's. Recasting 
this as a constraint on the seed and free solutions in terms of lower-order 
fields, we have 


n—1 

E 

1=0 


F„| + L„\h‘”’[k’l\,Fa\ 


t>0 




(61) 


Like the relaxed Einstein equation, Eq. (61) can be solved order by order 


in r, thereby constraining the functions I^^i^ and which in the general 
solution to the relaxed Einstein equation were entirely arbitrary. Working 
through sequential orders in r first constrains each in sequence, from 
£ = n — 1 to £ = {), and then each k"^^^ in sequence from £ = tlto £ = oo. The 
most important outcome is evolution equations for the moments I^yi^{t). 
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At n = 1, the term with the lowest power of r in Eq. (26a), l/r^ 

~ C>(l/r), which determines = 2m{t)6fj,,y. 
The overall factor is written as 2m{t) for good reason; from Eq. (|42|), the 

inner background metric has an expansion + 0 (l/r^), and sub¬ 

stituting that expansion into the formula for the Arnowitt-Deser-Misner 
(ADM) mass, one finds that m{t) is precisely the ADM mass of We 
call this the leading-order mass of the object. Proceeding to the next order 
in r, 1 /r, time derivatives and acceleration terms appear, and we find our 
sought evolution equations 


dm 

dt 


= 0 , = 0 , 


(62) 


which tell us that the object behaves approximately as a test body in g^^iy. 
The remainder of the content of the n = 1 gauge condition, at all orders 
in r, is to enforce various uninteresting relationships between the functions 


un 

At n = 2 , the lowest order in r in the gauge condition (26b), 1/r^, sim¬ 
ilarly determines (after a very slight gauge refinement that remains within 
the Lorenz gauge m 
the form 


that divides into mass and current moments of 


Ml,, = 2M,it)5, 




Slvi = 4«(^e^)ji5'^(t) 


(63) 


where euj = 0. M, and S'* are, respectively, the mass dipole moment and 
ADM angular momentum of g'^^] we call these the object’s leading-order 
mass dipole moment relative to 7 and leading-order spin about 7 . Defining 
5m^i/(t) = we also find (from the order-l/r^ and 1 /r terms in the 

gauge condition) that 


-h 4 m(o + 2Mp )), (64) 

where . 6mai3 can be thought of loosely as a correction to 

the object’s monopole moment, but it is not invariant; a smooth gauge 
perturbation that is nonvanishing on 7 will add directly to A:^^. This 

is why the moments are called algorithmic: they arise in the algorithm of 
solving the field equations, but they are not invariant descriptions of the 
perturbed object. Furthermore, is determined only up to a constant 

c5ai3, but I choose to incorporate that term into mdap- 

Finally, the 1/r^ term in the n = 2 gauge condition determines the 
evolution equation ^ = 0 , which tells us that the leading-order spin is 
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parallelly propagated along 7 , and the 1 /r term determines 

^^ 2 "" = -RaOboM’’ - mFl - + \RQaijS"\ (65) 

where I have dehned 5*-^ = I call this the master equation of motion. 

It describes how the small object moves relative to any choice of nearby 
worldline 7 . That motion is influenced by (i) geodesic deviation (evinced by 
the first term on the right), (ii) the acceleration of 7 (manifesting as T))), 
(iii) the ambient free helds in the neighbourhood, and (iv) the Mathisson- 
Papapetrou force created by coupling of the spin to the curvature of the 
external background metric. In the self-consistent expansion, we choose 7 
such that M* = 0. Therefore, 


dk^ 


= -- 


ta 


dt 


+ h^tta + ^^OaijS^R 


( 66 ) 


This is the MiSaTaQuWa force (albeit in an atypical form) plus the Mathisson- 
Papapetrou force. We see more clearly now that the practical consequence 
of the self-consistent approach, by utilizing a worldline relative to which M* 
vanishes, is to prevent the term from appearing in the second- 

order field (58). According to Eq. (65), this term would otherwise grow 


quadratically with time; this is the same undesirable growth discussed in 
See. 01 

As at n = 1, at n = 2 the gauge condition imposes relationships between 
the functions these relationships do not involve any of the algorithmic 
moments, and they leave each written in terms of arbitrary functions 
of time. 

Extrapolating to arbitrary order n, we draw several conclusions: (i) The 
moment appearing in the leading order term is directly 

identified with a moment of . Explicitly, as the mass appears at order 
and the spin at order e^, so the quadrupole moment will appear at order 

as 


e^, the octupole at order e^, etc ^ We interpret these moments 
the leading-order moments of the object, (ii) The others, 1”^,^ with i < 
n — 1 , we interpret as perturbations to the object’s moments, and the gauge 


^To relate back to the Thorne and Hartle results (|13l)-(|14l), canonical mass and current 


quadrupole moments Qij and Qij can be defined from MFfij and S'Ai, according to 


Mt 




S fjLuij j, 


(67) 


similarly to Eq. (63|. In choosing the normalization factors of these moments, I follow 

Ref. [70]. 
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condition determines that each of these corrections involve some linear or 
nonlinear factors of moments of g'j^^. (iii) Even after the gauge and matching 
conditions are imposed, every coefficient (t) appearing in the free fields 
is made up entirely of arbitrary functions of time that are determined only 
once global boundary conditions (BCl) are also imposed]^ 

When going to higher order, it is not obvious what condition should be 
imposed to make the worldline a good representation of the object’s bulk 
motion. A natural candidate would be to set the corrections to the 

object’s mass dipole moment to zero. However, one would have to investigate 
the third- and higher-order gauge conditions to see whether (i) doing so is 
possible and (ii) doing so prevents unwanted secular growth that can be 
associated with a growing displacement of the object away from the chosen 
representative worldline. Section offers a different method of choosing a 
good worldline, at least at second order. 


3.5 Split into self-field and effective field 

I now define a convenient split of into a self-field /i®” and an effectively 
external field Specifically, I define to be the piece of Eq. (59) that 
contains no linear or nonlinear combinations of the algorithmic moments 

jn . 
fii^L ■ 


h 


Rn 

flU 


Y, /ijir (^; Z, {kr}m<n,i>o), 

l>0 


( 68 ) 


where is the part of containing no factors of This implicitly 

defines the self-field /i®!l = /ijl,. — to be 


n—1 


hi: = E {ir}m<n,e<n, (69) 


e=o 


where every term in contains at least one factor of at least one al¬ 
gorithmic moment. I also define and 

proof of statement (iii) will be presented elsewhere. Here it can be taken as a 
conjecture, although it is known to be true at all orders in e and r that have been explicitly 
considered. Intuitively, it can be inferred from the fact that we can choose boundary 
conditions for which all vanish at a given n, and with that choice we must still be 
able to satisfy the gauge condition; hence, the constraints on the cannot involve 

the k'^^j^’s of the same n. 
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Conclusion (iii) at the end of the last section implies that each of these two 
fields separately satisfies the Lorenz gauge condition 


25 


More explicitly, at first order, with 


h% = z, m) = + O(r0), (70) 

k}) = + rkl^i{t)n^ + 0{r^). (71) 

£>Q 


Recall that for any particular solution of the form (47), the functions in the 
free helds are . 

At second order, with 


hfv = z, S) + /ip‘^(x; 2 ;, 5m) + z, m, {k}}e>o), (72) 

K7{x] z, kj) + z, {k\]^>o), (73) 

£>0 

where 

(a:; 2 :, S) = + 0{l/r), (74) 

h^-f{x-z,5m) = ^-^ + 0{r\ (75) 

h^l = k%{t) + rkl^i{t)n^ + O(r^), (76) 

{x-, z,m, {kj}e>o) is made up exclusively of terms quadratic in m or 
products of m and (and derivatives of h^^), and {x; z, {kj}£>(j) is 

made up of terms quadratic in (and derivatives of h^^). 

I remind the reader that in Sec. |1.1.2] l identihed three “nice” properties 
I wish an effectively external metric to satisfy: 


1 . it should be a vacuum solution 


2 . it should be causal on the worldline 

3. the object should move as a test body in it. 

There is no guarantee a priori that these conditions can simultaneously be 
satisfied at all orders in e, but with the definitions given above, the first nice 

the statements were not true, one could always slightly alter the singular-regular 
split to make the two fields independently satisfy the gauge condition. Doing so would 
involve appropriately moving part of the free fields into hi. m- 
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property is met: the effective metric is a smooth solution to 

the vacuum Einstein equation Rfj.u[9] = 0 even at r = 0. This can be seen 
from the facts that (i) by construction, for r > 0 we still have a solution 
even if we set all moments to zero (i.e., if there is no object in B^), and 
(ii) for the free fields the construction applies even at r = 0. Note that here 


I am taking a particular series (68) obtained outside the object and simply 


analytically extending it down to r = 0. In this way, I define an effective 
spacetime Ad) in B^, where the object lies in the physical spacetime. 

More explicitly, the first- and second-order regular fields satisfy the vac¬ 
uum equations 


77' 

-C//y 


L^[h^\Fo] = 0, (77) 

L^[h^^,Fo] = -5L^[h^\Fi], (78) 


for all r > 0. The self-fields are left to satisfy 

.[AS1]=0, L^[h^\Fo]=0, 




= 25^R^,[h^\ /^si] + 4S^R^,[h^\ 
L^[h^^Fo] = -5L^[h^\Fi]. 


(79) 

(80) 
(81) 


The third “nice” condition is also met, at least through order e: com¬ 
bining Eqs. © and (|66[), we get, in covariant form, 


D^zt^ 


= - Ka;u)u’’u- + ^R^pau'^SP- + O(e^), (82) 


,R 1 


2m 


where I have defined the antisymmetric 4-tensor = e^e^S''^ 


^ Together 


with ^ = 0, Eq. (82) is the equation of motion of a test body in It 
can be put in the form of the Thorne-Hartle equations (13)-(14) by using 
the results of Appendix A.l to absorb the terms into the covariant 
derivative and proper time in g^y, using = mu^ + 0{e^) and the fact that 


m is constant, and using Eq. (46) to write the Riemann tensor in terms of 
the tidal field Bij (noting 55^ = Bij + 0{e)). Doing so allows us to identify 
gfj_y with Thorne and Hartle’s “external” metric at least in the weak 
sense that gfj,y = -|- 0(e^, er^); since only first derivatives of h^y appear 


in Eq. (82), we can make no claims on agreement at order m or higher. 


“The monopole correction (641 can also be trivially rewritten in terms of as 
Smap = \m + rn{gap + 2uaUp)u'^u'h^l -|- , (83) 

where I have set M* = 0, and all fields are evaluated on 7 . 
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I show in the next section that my definition of the regular field also sat¬ 
isfies the remaining nice property: it is causal on the worldline. I also note 
that at least through order , my dehnition of agrees with the Detweiler- 
Whiting dehnition [Sjj, despite the fact that the Detweiler-Whiting held is 
dehned in a wholly different manner based on Green’s functions PES]. How¬ 
ever, despite all these reassuring facts, I stress that in general, even if the 
effective metric I define is found to satisfy each of the three “nice” proper¬ 
ties at a given order e^, it is not the unique field satisfying those properties. 
One can simply shift any given free held z,kf ) (with n' < n and ^ 

sufficiently large for appear in the nth equation of motion), and 

its nonlinear combinations, from the regular held into the self-held. The 
strongest of my “nice conditions” appears to be that the object behaves as 
a test body in the effective metric: since higher multipole moments couple 
to higher derivatives of the effective metric, this strongly constrains which 
part of the full metric should go into the effective one. At any hnite per¬ 
turbative order, however, one can always alter the effective metric’s higher 
derivatives without spoiling any of its nice properties. Nonetheless, I deem 
my “nice” choice most natural as a part of the process of solving the relaxed 


Einstein equations using the local expansion (47): before making reference 


to global boundary conditions, I simply put all the terms that involve the 
object’s multipole moments into the self-field, and I put all the terms made 
up entirely of unknown functions into the effective field. 


3.6 Summary at first and second order 

To extract the key pragmatic information of the preceding sections, I restate 
the conclusions at first and second order: The first-order held is given by 
Eqs. where the self-held /i®), is locally determined by the (con¬ 

stant) mass m, and the effective held h^^ is to be determined by global 
boundary conditions. The second-order held is given by Eqs. (72)-(73), 


where the self-held h^ is locally determined by (i) the mass m, (ii) the 
hrst-order effective held (hi) the (parallely propagated) spin and 
(iv) the monopole correction given in Eq. (83); the effective held h^^ 


is to be determined by global boundary conditions and by h^^. I encourage 
readers to examine the more explicit expressions for these helds given in, 
e.g.. Ref. 


Finally, the object moves on a worldline governed by Eq. (82). 
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4 Algorithm for an nth-order self-consistent ap¬ 
proximation: point particles, punctures, and global 
solutions 


With the general solution in the buffer region ready at hand, I now describe 
how to use it to generate a global solution. In short, this relies on a rigorous 
procedure of replacing the physical field in B.y with a fictitious field without 
altering the field elsewhere. What are the fictitious fields? Simply the ana¬ 
lytical expansion (59) from the buffer region, extended to apply to all r > 0. 
This was done already for h^^^ in the previous section. Continuing the ex¬ 
pression for into makes it into a field that diverges at r = 0 : the true 
self-field in the interior of the body, whatever it may be, is replaced with 
this divergent field, and the self-field becomes the singular field. Since this 
extension into B.y does not affect the field in the buffer region, it also does 
not affect the field values outside the buffer, out in the external universe. 

Section 4.1 describes how this analytical extension can be used to ascribe 
a certain pointlike stress-energy distribution to the object, thereby recov¬ 
ering the point particle picture at linear order (though this stress-energy I 
construct cannot be taken as a physical source in the Einstein equations be¬ 
yond linear order). Section 4.2 then describes how the analytical extension 
can be used to obtain a global solution at any order; this will be the first 
point at which a global boundary condition (BCl) is finally imposed. 


4.1 Skeletal stress-energy 


Consider the solution (59) extended down to r > 0. For reasons described 


in Sec. 1 . 2 , there is no known distributional source for this solution on a 


domain that includes r = 0. However, we can devise the following setup: 
each of the seed solutions can be thought of as being sourced by a pointlike 
stress-energy supported on 7 . Everything else in the general solution then 
grows from these seeds, either being generated by nonlinearities or in the 
case of the free fields, being determined by global boundary conditions. 

I refer to the stress-energy for the seed solutions as the skeletal stress 
energy, taking after a similar phrase in Ref. m-, the idea is that in some 
sense, the object (or body) can be replaced by a skeleton. That skeleton is 
made up of multipole moments living on the object’s worldline. 

For each seed solution z, If), I roughly follow the approach taken 

by Gralla and Wald at first order [20], defining the (trace-reversed) distri- 
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butional stress-energy tensor 


T;t = ( 84 ) 

The right-hand side can be written as 

jn 

(85) 

where N= d^di 

First examine the most singular term. Using the identities dLr~^ = 
{—lY{2i — 1 )!!^^ and d'^dir~^ = — 47 rh^(x), we have 


TT 

d^di 


n 




4vr(-l)^+il-, 

-(9l( 5 (x). 


( 2 £- 1 )!! 

Integrating the right-hand side against a test function we find 


V'"I'^vLdLS^{x)yf^d^xdt = {-if / di dt 


= (-1)' / 




dt^ 


( 86 ) 

(87) 

( 88 ) 


where in going from the second line to the third I have utilized the identity 
= ^^iV~9- Here I have defined tensor that 

agrees with if all a* are spatial indices and zero otherwise, meaning 
that -g grpp respect to g^u and that = 0 . 

Equation ( 88 ) shows that 


I[f)dLd^{x)= f I^ff^"'‘^^S{x,z).a:,-aidT, 

J-i 


Equation (85) now reads 


E, 


/ii/ 




(89) 


(90) 


Now note that, by construction, vanishes pointwise for r > 0. 

If it is nonvanishing as a distribution!^ it must have support only on 7 , in 


'That it must be well defined as a distribution follows from it being the result of linear 
operations on and itself being a sum of terms constructed from linear 

operations on an integrable function. The latter fact follows from the first term in the sum 
being expressible as the linear operation Ol on an integrable function proportional to r~^ 
[as in the text above Eq. ( 861 ], and all higher-order terms in the sum being constructed 
from linear operations on lower order terms in the sum (as described in Sec. 3.2 1 . 
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which case it must be proportional to or a derivative thereof; but from 
the calculation just performed, that would lead to a different algorithmic 
moment (with ^ t) appearing in the solution, which would contradict 
the definition of the seed field So vanishes as a 

distribution. 


Therefore, Eq. (84) becomes 

(-1)' 


rr(x;z,/,") = 




(91) 


4(2£- 1)!! 

If we add up all the multipole moments, we arrive at a skeletal stress-energy 
tensor 


n,l I ■ 




(92) 


where I have defined the corrected and normalized moments = 

Sn ■ Notably, this skeletal stress-energy agrees in form with 

that of the traditional multipolar expansion of a material body’s stress- 
energy tensor [371 [mils]. 

At first order in e, there is only one seed field, z,m), and Eq. (91) 

gives 

(93) 


T(^^(x;z,m)= / mu^u‘'6{x, z)dT 

J-y 


1.1 


Because = 


—precisely the point particle stress-energy used in Sec. 

x;z,m), and = 0, we can conclude that = 

T^ixiz). In other words, the first-order field is identical to one sourced by 
a point particle. Of course, this is derived from the analytically extended 
held, but it also applies to the physical held at distances r S> e from the 
worldline. 

At second order, there are seed helds, z, 6m), h^^^‘^{x; z, S), and 

z,M)-, although the last of these three we set to zero with our choice 
of z^, it is worth displaying the skeletal stress-energy that would source it 
if we chose a different z^. Equation (91) gives 


1 


{x\ z, 6m) = T / dm^ 6{x,z)dT, 

^ J'‘y 

T^^{x-,z,S) = - f u^^S'^^^Vc,6{x,z)dT, 
J-f 

T^’'{x;z,M) = - [ u>^u'^M'^Vadix,z)dT, 


(94) 

(95) 

(96) 
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where the overline indicates trace reversal. Appendix |B.l| shows that zq, M) 

[plus a piece of T!^'^[x] zo,5m)] is precisely equal to the linear term in the 
expansion of z, m) given z^ = Zq + eM^/m + 0 (e). 


4.2 Puncture schemes 


I now describe how one can use the results from the buffer region to obtain 
a global solution. The method is called a puncture scheme. It has become 
standard in the linearized problem |56H59] . and it is the only known practical 
way to obtain numerical results at second order and higher [m [m Ea [sa 

[73]. 

A puncture scheme begins with the construction of a puncture de¬ 
fined by truncating the local expansion of the singular field from the buffer 
region at some specified order. One then defines the residual field 


= h — 

"'IMU — "'IlL/ 


(97) 


and in a region covering the object, writes a field equation for rather 

than one for (the analytically continued) physical field Since « /i®^, 

so too The better represents the better represents 

For example, if lim^: 


- ^L(^)] = 0 > then = 


"fiv ^XX xxxxxa;_>xy|^,t^j^i^.^y fj,uJ1 ~ 'y ‘ ^^1/ 

lima;_s.xy h^^{x)\ that is, the residual field agrees with the regular field on the 
worldline. If agrees with to one order higher, meaning = 

o(r), then = \\m.x^-yVph^^] since the self-force is constructed 

from first derivatives of this condition guarantees that the force can be 
calculated from as in Eq. (100) below. 

There are several schemes that can be developed from the starting point 
of the puncture. Here I describe a worldtube scheme in the tradition of 
Refs. [SniESl- In this type of scheme one uses the field variables /ijf” inside 
a worldtube T surrounding 7 , the field variables outside that worldtube, 
and the change of variables when moving between the two 

regions,^® Concretely, a second-order puncture scheme is then summarized 


^®One does not solve the problem in each domain separately, since the separate problems 
would be ill-posed. Instead, when calculating at a point just outside F that depends on 
points on past time slices inside F, one makes use of the values of already calculated 
at those earlier points, and vice versa; see Sec. VB of Ref. |56 |. 
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by the coupled system of equationj^ 
Ei,y[h^] = 0 


inside F, 
outside r, 


(98a) 

(98b) 


E^,[h^^] = 25‘^R^,[h\h^] - E^,[h^^] = S\ 
E^,[h^]=25‘^R^,[h\h^] 


eS2 

fiV 


inside F, (99a) 

outside F, (99b) 


^ is.-’ - /.?■') - ftSs„) 


u u 


/9 


( 100 ) 


where the puncture diverges on the worldline determined by Eq. (100). 


Here I have looked ahead by using the second-order equation of motion 
(for an object whose spin and quadrupole moments vanish); the first-order 
equation could be used instead, though as discussed in Sec. |1.5.4[ the 
results would be accurate in a smaller region of spacetime. 

-( |100[ ) must be solved together, as a coupled 
system for the variables (inside/outside F), and The 


In this scheme, Eqs. (98) 


residual fields govern the position of the puncture, and the position of the 
puncture effectively sources the residual fields. This system of equations is 
to be solved subject to some global boundary conditions. For simplicity, we 
can consider using specified initial data on a Cauchy surface 


I make two techical asides: 


according to Eq. 


72 


contains terms 


proportional to the first-order regular field So at each timestep in a 

numerical evolution, one must first calculate the first-order residual field 
from the first-order puncture, then use that residual field to calculate the 
second-order effective source, and so forth. Also, the divergence of the indi¬ 
vidual sources is quite strong at the worldline, with the terms h^] 

and in Eq. (99) each going as l/r"^; by construction, these diver¬ 

gences analytically cancel each other, but the cancellation is numerically 
delicate. 

In principle, there is no obstacle to using a puncture scheme at any 
order in e. Outside F, one may solve the physical problem Efj,,y[h^] = 


^®The effective sources are usually written to include the skeletal stress-energy 

terms, which are canceled by distributional content in Here I have instead 

followed Ref. | 22| by writing the source pointwise, off 7 ; if the puncture agrees with the 
singular held sufficiently well, the source at points on 7 can be dehned as the limit from 
off 7 . 

Because the approximation is accurate only within a hnite region of size 1 /e, one 
might better solve the equations in the future domain of dependence of a partial Cauchy 
surface. 
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S^^[h ^,..., inside F one may solve the effective problem = 

S^^[h ^,..., — Ef^,y[h'^'^] = 5'®®”', and when crossing F, one may change 

variables from the residual field to the full field via + /i^”. Of 

course, this requires sufficiently high-order expressions for the puncture and 
for the equation of motion governing how the puncture moves. But we know 
how to obtain both expressions from a local analysis in the buffer region. 

The basic picture of a puncture scheme is illustrated in Fig.[^ Physically, 
the scheme replaces the actual problem in the region covering the body, 
with all its matter fields, singularities (in the case of a black hole), or other 
oddities (in the case of exotic matter), with an effective problem. But it 
yields the correct physical field outside the object. Hence, while we must 
abandon the point particle model at nonlinear orders, we can replace it with 
the more general concept of a puncture, a local singularity that encodes all 
the necessary information about the object in its multipole moments. 


4.3 The causality of the regular field 


I now return to the question of how “nice” the regular field is, by which I 
specifically mean the three enumerated properties in Sec. 3.5 I have already 
shown that at all orders, the regular field I defined is a vacuum solution. I 
have also shown that at first order, the object moves as a test body in the 
effective geometry it induces, and (for spherical, nonspinning objects) this 
result will be extended to second order in Sec. Now the final property 
follows immediately, at all orders in e, from the design of the puncture 
scheme: the residual field and its derivatives manifestly depend only on the 
causal past, and by construction, the residual field and any number of its 
derivatives agree with those of the regular field on the worldline. Taking the 
limit of an infinite number of derivatives, we get the desired result. 


4.4 A note on “regularization” 

In the gravitational self-force literature, one often speaks of “regularizing” 
the field or the self-force. This can mistakenly give the impression that 
one has introduced infinities into the problem, and that one must regularize 
them to recover the physical result. But in the formalism I have described 
here, one only ever deals with finite quantities. The force on the object, 
for example, is derived from the field equations outside the object, and it 
is written in terms of manifestly finite fields in that region. Those fields 
outside the object are then written in terms of quantities on the worldline 
by identifying 5*(f), a field off the worldline, with e^^S^{z{t)), for example; 
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or similarly, by defining the regular field such that h^^{z{t)) = and 

From this perspective, the various “regularization” methods that have 
been used in the self-force problem to remove the ‘singular part’ of the 
field [8l |9] arise only as a practical necessity: we cannot easily determine the 
physical metric inside the object, nor are we interested in doing so, which 
prompts us to replace it with the fiction of a singular field solely as a means 
of calculating the physical metric outside the object. Computational tech¬ 
niques such as puncture schemes and mode-sum “regularization” BM are 
not methods of removing singularities; they are simply methods of calculat¬ 
ing the particular finite quantities in question. In mode-sum regularization, 
for example, one rigorously writes a spherical-harmonic mode decomposition 
of h^^{z) by decomposing h^l{z) = with being the 

held of a point mass. Every quantity in the calculation is hnite every step 
of the way. 


Gralla-Wald and osculating-geodesics approxi¬ 
mations 


Given results in the self-consistent approximation, one can always obtain 
analogous results in a Gralla-Wald or osculating-geodesics approximation 


by performing an expansion of the worldline, as described in Secs. 1.5.5 and 


1.5.6 In this section I present that procedure and those results. 


5.1 Gralla-Wald approximation 

There are two ways to obtain the Gralla-Wald approximation from prior 
results: by expanding the worldline in the self-consistent results; or more 
simply, by returning to the results in the buffer region and simply choosing 
the worldline at r = 0 to be 70, the zeroth-order worldline. 

For the moment, let us take the second route. All the local results are 
unchanged, except that mass dipole moment M* 

is not set to zero. This means that the metric perturbations 


are modihed from Eqs. (70)-(73) to be 


= h: 


seed 


{x-,zo,m) = 


2m6, 


fliy 


+ O(r0), 


= '^h^/lt%x]Zo,kj), 

e>o 


( 101 ) 

( 102 ) 
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full metric g 


fll/ 


"self field" 


/i® 

"'Hiy 


effective metric 





Figure 3: Replacement of the physical metric with an effective metric plus a 
puncture. In the top row, the physical metric is split into a self-field 
plus an “effectively external” metric In the bottom row, 

we “black out” the region in and very near the object, and we replace the 
physical metric with a singular field (any local approximation to which 
is called a puncture) plus the “effectively external” metric. The replacement 
of the self-field with the singular field is made only very near the object: the 
self-field and singular field are identical to one another in the buffer region 
and beyond; the effective metric g^y is completely unaltered everywhere it 
is defined. 


and 


zo, M) + zq, S) + h^^y'^{x] zq, 5m) 

+ (103) 

^0’ ^0, (104) 

£>0 

where 

O l\/f • Y) ^ 

zo, M) = -^5^y + 0(l/r), (105) 


and 5m^y differs from Eq. (83) by the inclusion of the Mi term in Eq. (64), 
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becoming 


5map = \rn { 2 h^l + + m{ga 0 + 2 uaUp)u^"v!'h^l 

+ Aui^o,{mh^)^u^^ + 2Mfi)). (106) 


Otherwise, all functionals remain completely unchanged, except that they 
are evaluated as functionals of Zq rather than of z^. The regular field nat¬ 
urally remains a vacuum solution even on 70 , its first two orders satisfying 


= 0 , and L^[h^^] = 0 . 

Finally, if we define z^_^ = M^/m to be the deviation perpendicular 
to 7 o (defining M* = 0), then from the master equation (65), we get the 
Gralla-Wald equation 


dT^ 




(107) 


Just as in the self-consistent case, these local results can be used in 
the design of a puncture scheme to obtain global results. However, before 
describing that scheme, I return to the other method of obtaining the Gralla- 
Wald approximation: by substituting an expansion of the worldline into the 
self-consistent results. 


5.1.1 Expansion of the worldline 

We would like to express the expansion covariantly, in terms of vectors 
that live on 79 . Suppose we begin with an expansion z^{s,e) = Zq(s) + 
€z^(s) + e^z^(s) + ..where z{i(s) = linear term, z^(s), 

is automatically a vector on 70 , since it is a first derivative along the curve 
of constant s and increasing e. However, beyond linear order, the terms 
are no longer tensorial; each of them is simply a collections of four scalars 
dependent upon the particular chart in which z^(s,e) = x^('y(s, e)). 

So let us approach the problem more geometrically. To facilitate the 
expansion, I introduce a Lie derivative £ that acts on a functional’s depen¬ 
dence on z^: 

£^Ap^{x]z) =-^AK{x-,z + \i)\^^^, (108) 

where A is a multi-index of any covariant and contravariant rank. This is 
closely related to a Lie derivative acting at £ moves the field point 
relative to the worldline; £ moves the worldline relative to the field point. 
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The two operations are not identical, since the tensorial character of the 
functional is different at the two points. For example, in Eq. (217) we see 
that z) is a rank-two tensor at but a scalar at z^. 

An expansion of a functional of the worldline z^{s, e) in the limit e —)• 0 
is really an expansion along the flow of increasing e. We can write this as 

A^{x\ Ze) = Aa(x; zq) + e5A^{x-, zq) -h e^(5^AA(x; zq) + • • •, (109) 


where 


(5’^Aa(x;zo) = -7"^^”-4a(x;zo); 


ni 


( 110 ) 


here is the generator of the flow. 


Now, let Z 2 p = This is a vector on jq. It is equal to z^ 

if the expansion of the worldline is performed in Fermi normal coordinates. 
One can obtain a covariant evolution equation for it from a second-order 
self-consistent equation of motion, as described in Appendix |A.2[ Beginning 
from Eq. (191), this procedure leads to Eq. ( 205[ ), for example. More im¬ 
portantly for the present analysis, we can also write the expansion (110) in 


terms of this quantity. As an instance of that, the second-order term is 


]^£lAt^{x] Zq) = z^^V^/Aa(x;zo) -h Zi'v^/Vi./Aa(x; zq). 


( 111 ) 


where primed indices refer to the point Zq . One can do the same at any 
order: from the equation of motion for z^, obtain evolution equations for 


vectors on Zq, and write the expansion (109) in terms of those vectors. 


Applying this expansion to the perturbations from the self-consistent 
approximation, we have 


e^hl,{x;z) = e^hl,{x;zo)+e^+^£,,hl,{x-zo) + 0{e^+^). ( 112 ) 

Therefore, the Gralla-Wald expansion ( [^ reads 

= eh]^u{x-, Zq) + e^h^^ix; zq, zi) 0{e^), (113) 

with hjj^j^{x] Zq) being the same functional as hj^^{x;z), but with Zg having 
replaced z^ in its argument, and with the second-order perturbation given 
by the new functional 


Zq, Zi) = hl^{x] Zq) + Sh^^ix; Zo, Zi). (114) 

where ShJ^^(x; zq, zi) = £z^hf,^(x; zq). 
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The same expansions are applied in precise analogy for the singular and 
regular fields, yielding with 


= h%^^{x-zo), (115) 

zo) + 6h%^^{x- zq, zi). (116) 


Appendix describes these expansions in more detail. But in a few 
words, the end result is that the procedure exactly recovers Eqs. (101)-(107). 
All of the meat of the result is in 


{x; zq, zi): the action of £ on h] 


tiu pre¬ 
cisely generates the mass dipole moment seed field zq, M) and the 

contribution of the mass dipole moment to 5m Appendix also shows 
that (i) the generation and modification of these two seed fields precisely 
corresponds to the linear term in the expansion of the skeletal stress-energy 
z) around z^ = Zq , and (ii) the term 6hjj^^{x-, zq, zi) [and analogously, 

5y?J,^^{x\ Zq, Zi)] can be written as a functional 5h}^^{x\ zq, zi±) that depends 
only on the perpendicular piece of z^. Result (ii) is in agreement with the 
fact that the seed fields depend only on = mz^j_. Physically, this results 
from the field equations’ indifference to the piece of the deviation that is tan¬ 
gential to the worldline, which can always be set to zero by reparametrizing 
the family z^{s, e) with a parameter s' = s'(s, e) that ensures curves of fixed 
s' strike Zq orthogonally at e = 0 . 


5.1.2 Global metric 

A global solution can be obtained from the local results just as in the self- 
consistent case. I first note that the local singularity structure of the singular 
field = ehfl -k is identical to the self-consistent singular field, but 

for two important alterations: 

• The divergent terms diverge on 70 , not on 7 . 

• The second-order singular field depends on the correction z^ to the 
position. 


Because the point at which the puncture diverges is independent of the 
perturbations in this expansion, the puncture scheme becomes a se¬ 
quence of equations, rather than a coupled system: first, the zeroth-order 
worldline is prescribed as a solution to the background geodesic equation. 




= 0 , 


(117) 
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then the first order field is found from 




PI] 


Ef,u[h^] = 0 


inside Fq, 
outside Fq, 


(118a) 

(118b) 


then that field is used to find the first-order correction to the position by 


solving the Gralla-Wald equation (107), and finally the second-order field is 
found from 


= 2S^Rf,,[h\ h^] - E^4hP^] 
Ef,4h^] = 26^R^^[h\V] 


inside Fq, 
outside Fq. 


(119a) 

(119b) 


Here Fq is a tube around 70 ; unlike in the self-consistent case, neither 70 nor 
Fq need be updated over the course of the numerical simulation. Like the 
self-consistent puncture scheme, these equations can be solved given data 
on a Cauchy surface 
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Again, there is no obstacle to carrying a puncture scheme like this to 
arbitrary order, given the local solution in the buffer region. 

5.2 Osculating-geodesics approximation 


I refer the reader back to Sec. |ig for a reminder of how the osculating- 
geodesics approximation works. In brief, it obtains a self-consistent ap¬ 
proximation by first using a sequence of Gralla-Wald approximations to 
find the the self-consistent worldline z^, then solving the relaxed field equa¬ 
tions with already determined. Concretely, through second order one 


seeks a solution to the self-consistent system (98)-(100) by applying the 


Gralla-Wald approximation (113) to the right-hand side of Eq. (100) at 


each instant r, with the expansion z^{T',e) = <(r)('^') + + 0{e) be¬ 

ing around the geodesic that is instantaneously tangential to z^ at 


time r. The terms E^^{z{t); zq, ..., Zn-i) that then appear in Eq. (100) are 


found by solving Eqs. (117)-(119). By using the resulting sequence of forces 


eF(^(r; Zo(t)) Ee'^E^ir; 2 ;i(t-)), one can solve Eq. (100). Einally, one can 
solve Eqs. (l98|)-(99) using the z^ one has already found. 


What is the advantage of this over simply solving Eqs. (98)“(100) di¬ 
rectly? At linear order, the advantage is that one can solve Eq. (|118) 


very easily in the frequency domain, while the self-consistent equations 


Since the approximation is held to be valid in a region of size cfe) ^ 1/v7, a reasonable 
approach would be to solve the equations in the causal future of a partial Cauchy surface 
of that size. 
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do not seem to admit a frequency decomposition. One can then com¬ 
pute a table of values of (r; 2;o(r)) different geodesics and easily solve 
•^^2^ = eF|‘(r; Zo(t)) to obtain an approximation to [62]. But at second- 
order, it is not obvious whether the scheme is simpler than a direct solution 
of the self-consistent equations. The term Zo(t-), •?i(r)) that appears 

in the second-order field (and a piece of which appears in the puncture h^) 
grows with time away from the osculation point can it be han¬ 

dled in an efficient way? Furthermore, one must ensure that appropriate 
boundary conditions are used in each Gralla-Wald expansion to reproduce 
desired global boundary conditions (such as no incoming radiation) on the 
self-consistent field; since the approximations only apply in a region of size 
<C l/\/e, is there any simple way to find these boundary conditions? I leave 
these questions open. 


6 Gauge transformations in perturbative descrip¬ 
tions of motion 


Thus far, all the explicit results I have presented have been confined to a 
single choice of gauge. But as I described in Sec. |1.4[ gauge and motion 


are intimately related in perturbation theory: any worldline z^ one finds 
can always be shifted by an amount Precisely how this impacts one’s 
approximation scheme depends strongly on how one represents the worldline. 

The effect of a gauge transformation on the worldline (or equivalently, on 
the self-force) was hrst explored by Barack and Ori |^. Their results were 
extended to discontinous gauge transformations by Gralla and Wald [201 
I2lj^ and to even more singular transformations in Ref. |33j . In this section 
I focus, for simplicity, on smooth transformations, and I aim mostly at 
(i) clarifying the gauge freedom in each of the different representations of 
motion and their corresponding approximation schemes, and (ii) dehning 
appropriate rules for the gauge transformations of the singular and regular 
fields. The presentation is adapted from Ref. |19j . which includes additional 
results and proves some key statements that are here asserted without proof. 


^^The extension of Gralla and Wald’s result to the self-consistent case EU unfortunately 
contained a significant error, leading to a result that held only in gauges continuously 
related to Lorenz, as in the Barack-Ori analysis. 
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6.1 Review of gauge freedom in perturbation theory 

Before considering the question of gauge in self-force theory, I briefly remind 
the reader of the basics, following Refs. [71H76]. It will be convenient in this 
section to adopt index-free notation for tensors. In perturbation theory, we 
consider a family of metrics g^iy(x, e), or simply g in the absence of a chart. 
This family lives on a family of manifolds Ate, and a given choice of gauge 
refers to an identification map cj)^ : Ato Ale- The identihcation map 
induces a flow through the family down to the base manifold Ato where the 
background metric g = go lives. Call the generator of this flow X = 

We wish to approximate a tensor A at a point (t>f{p) G Ate as an expansion 
around its value at e = 0. This expansion is given by 

fTl 

iXm) = (e"'-‘A)(p) = ^ -y iCkm), (120) 

n>0 

where (j)^* denotes the pullback of , C is the Lie derivative, and p G Mq. 
We define the nth-order perturbation in this gauge to be 

Aj{p)^i(£lA)(p). (121) 

n\ 

Now say we work in a different gauge. This corresponds to a different 
choice of identihcation map : Ato —>■ At^ and how generator Y = 

The approximation of the tensor A in terms of tensors at the point p G Afo 
is now given by 

{4>^*A){p) = {e^‘^^k){p), ( 122 ) 

and the nth-order perturbation is 

(123) 

n! 

The gauge transformation of A is the difference between the two expan¬ 
sions when evaluated at a point in AIq^ in more common notation, we say 
An —> A'n = An + XAn, where the primed tensor refers to the Y gauge, the 
unprimed to the X gauge, and their difference is 

AA„(p) = G(2;fA)(p) - G(£’iA)(p). (124) 

n! n! 

The hrst- and second-order terms are easily expressed in the familiar form 
AAi = C^^Ao, (125a) 

AA 2 = C^^Aq + —Ai, (125b) 


62 


where = Y — X and ^2 = usual gauge vectors. Higher- 

order terms are straightforwardly written down, but since concrete self-force 
results are not available beyond second order, I stop here. 

In a chart, one can show that a gauge transformation can equivalently 
be thought of as a near-identity coordinate transformation rather than a 
change in identification map. First lay a chart on each using some 
identification, say {p)) = x^{p) for each p G A^o- The two identifica¬ 

tion maps (j)^ and identify the point p in Ado with two different points 
q = (p) and q' = cj)^ {p) in Ade, which are related by the active diffeo- 

morphism q 1 —)■ (pj{{4'^)~^{q)) and which have slightly different coordinate 
values. Since the chart x^ consists of four ordinary scalar fields, we can 


apply the general expansion (122) to write the coordinates at q’ as 


x^{q') = x^{q) + eCf(x(g)) -h ^2 + \^iix{q))d„^^{x{q)) 


+ 0{e^ 


(126) 

where I have used Cxx^ = 0 to express Y derivatives as ^ derivatives, 
and I have used the fact that x^{p) = x^{q) to express the components 
on the right-hand side as functions of x^{q). Now, rather than an active 
diffeomorphism on Ad^, let us consider this as a passive change in the chart, 
x^{q) I— x'^{x{q)). Define the coordinate transformation such that x'^{q') = 


x^{q). Rewriting Eq. (126) as an equation for x^(g) as a function of x^(g'), 
we get 


x'^{q) = x^^{q') - (x(g')) - ^ 2 (^( 9 )) “ \ii{x{q))duii{x{q)) 


+ 0(e3). 


(127) 


Gauge transformation laws for components of tensors can be derived 
directly from this coordinate transformation. For example, by rewriting 
Eq. (127) as an equation for x^{x'{q')) and substituting it into the ordinary 
transformation law for the components of the metric, one finds 


3t^ 

= g/iAx', e) + Y^gf,u{x', e) + e) O(e^), 


(128a) 

(128b) 


where the analogous transformation law for a 

tensor of arbitrary rank is also easily found. To relate this to the language 
used above, note that we are now using a single identification map and 
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given that identification, components of g have identical coordinate values 


as {(pe*g)- Equation (128) applies even if g^y{x', e) is not expanded for small 


e. 


If it is so expanded, then Eq. (128) returns Eq. (125). 


Of course, we are ultimately interested in solving the Einstein equation, 
and so in addition to transformations of the metric, we must consider trans¬ 
formations of curvature tensors. Useful identities for the transformations 
of curvature tensors are derived in Appendix One of their consequences 
is that when examining perturbations of a curvature tensor, one can derive 
transformation laws in two equally natural ways: directly from Eq. (125) 


or from the transformations of the metric perturbations. For example, in a 


vacuum background, Eq. (125) directly implies 


AdR^u[h^] = C^^R^^lg] = 0 , 

A{5Rf,,[h^] + 5^R^,u[h\ h^]) = C^6R^,[h^]- 


(129) 

(130) 


or the same equations can be found by instead using Eq. (125) for the metric 
itself, writing 


A{6R^,[h^] + 6^Rf,,[h\h^]) 

= 6R^,[h'^] + 5^Rf,,W\ - {SR^Ah^] + S^R^,[h\h^]) 

= 5R^u[^^ 29] + + SRfj,u[C^^h^] 

+ + 6‘^R^|2[C^-^^g, C^^g] (131) 


and then applying Eqs. (231)-(233). 


6.2 Gauge in the Gralla-Wald approximation 

I reverse my usual ordering by first considering transformations in the Gralla- 
Wald approximation, which allows the most straightforward treatment. Since 
the Gralla-Wald expansion is an ordinary one, with coefficients independent 
of e, all the ordinary rules apply. 


6.2.1 Transformation of the metric and the worldline 

First, let us examine the transformation of the deviation terms in the expan¬ 
sion of the worldline. According to Eq. (127), under a gauge transformation 
the coordinates z>^{s,e) = x^(%(s)) on the worldline become 


z'^(s,e) = z^(s,e) -eCf(2) ^^(z) - ^^i(z)d^^^(z) 


+ 0(e^), (132) 
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where functions of are evaluated at z^{s). By expanding this in powers 
of e, we immediately find 


(133a) 

z[^(s) = z^(s) (133b) 

4^(s) = 4(4 - 4(^o) + ^4(zo)du4(^o) - z'({s)du4(zo), (133c) 


where functions of Zq are evaluated at Zq(s). Note that the zeroth-order 
worldline is unchanged; the effect of the transformation is to alter the devi¬ 
ations relative to that worldline. 

Now let us examine the transformation of the metric perturbations. Ap¬ 


plying the transformation laws (125) to the metric, we find 
with 






(134a) 

(134b) 


Corresponding to the treatment of the worldline, all like all /i”^, diverge 
at Zq(s). Instead of altering the curve on which the fields diverge, the gauge 
transformation alters the singularity on that curve, by altering the functions 
that appear in 


In fact, one need not even appeal to Eq. (132) to determine how the de¬ 
viation vectors are transformed: the transformation laws for z((>g, as given 
in Eq. (133), can be derived directly from those for For example, z^ 


appears in the metric as the mass dipole moment M^/m, and the trans¬ 
formation of can be determined from the transformation law for 


By applying Eq. (134b) to the results (101)”(104) in Fermi normal coordi 
nate: 


and comparing to Eq. (105), one quickly finds that is the 


only part of that contributes a mass dipole term. More specifically. 


2m5, 


_27T2-^^ 71 / ‘ 

— contributes the relevant term, given by-This modifies 


%■ 

IvP by an amount AM* = —m^\, from which we read off the same result as 


in Eq. (133b). 


6.2.2 Transformation of the singular and regular fields 


Though deriving the transformation laws for was trivial, we must put 
some thought into how to do the same for the singular and regular fields h®” 


this stage the transformation is applied only to the sums + h®™; 

does not yet need transformation laws for the individual pieces and 
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and We are free to split into a singular piece and a regular piece 
however we like, just as we were with However, ideally, we can do so 
in such a way that we preserve the nice properties of the split. Allow me 
to specialize to a nonspinning, spherical object, such that there is no force 
from coupling of moments to external curvature. This means the properties 
I wish to preserve are that the transformed worldline z'^ [in its expanded 


form (133)] should satisfy the geodesic equation in the transformed effective 
metric g'^^, = and that effective metric should satisfy the vacuum 

Einstein equation. 

Appropriate transformation laws can be found by noting that for a 
smooth metric and smooth worldline, both the geodesic equation and the 
vacuum Einstein equation are manifestly invariant under a generic smooth 
coordinate transformation. It follows that if the worldline and the metric in 
those two equations are expanded in a power series, and the metric and the 
expanded geodesic transform according to the standard laws of gauge trans¬ 
formations, then the perturbative expansions of the worldline and metric 
continue to satisfy the perturbative geodesic equation and the perturbative 
Einstein equation. In our case, we already have the worldline transform¬ 
ing according to the standard transformation law. Accordingly, if we de¬ 
fine the regular field to transform like an ordinary smooth perturbation, as 
= + where 


/^h^},{x]zo) = C^^g^u, 


Ah 


R2 


(x, .^o) T T (x, .sq)? 


(135a) 

(135b) 


then all the nice properties of the regular field are maintained under a gauge 
transformation. However, at the same time, the total field h^y = -(- 


must satisfy Eq. (134). This means the singular field must satisfy 


-or more explicitly. 


= hi" + - A/i 




Rn\ 

fii/ ) 


Ahll{x-,zo) = 0 , 

Ahll{x;zo) = C^^hllix] zo) 


(136) 


(137a) 

(137b) 


6.2.3 Governing equations in alternative gauges 


By design, the transformation laws (133)~(137) ensure that the governing 


equations of the Gralla-Wald representation are invariant under a gauge 
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transformation: in any gauge, the regular field satisfies the vacuum equa¬ 
tions = 0 , SRij,,y[h^^] = Zq is unaltered and 


the deviations from it satisfy (|35|)-(|36|):^^Bnd assuming the gauge choice 


admits a well-posed initial value problem, the full field off the worldline 


can be calculated using the puncture scheme encapsulated in Eqs. (117)- 


(119). These facts follow from the invariance of the geodesic equation and 


Eqs. (129)-(131). However, because of Eq. (137b), the singularity structure 


of the second-order field, and hence the puncture in the puncture scheme, is 
altered by a gauge transformation. We can assess the impact of this change 
by examining how it manifests in the second-order field equation. Using 


Eq. (233), we have 




Sli 


,S1 




(139a) 

(139b) 


The second term on the right is precisely the term required to make a so¬ 
lution to the correct equation, 

at points off 70 . The first term on the right-hand side of Eq. ( |139b ) indicates 
a change in the skeletal stress-energy: —)• This change 

in skeleton corresponds to a change in the seed solutions ZQ,6m) 

and zo, M). An explicit expression for can be derived from 

one for , given in Eq. (215). From that latter equation, we see that 

both alters 5m^j,y and more notably, shifts the mass dipole moment 
by an amount AM^ = —in agreement with the discussion above. 

Of course, to make use of the fact that the governing equations are the 
same in all smoothly related gauges, one must first find an effective metric 
satisfying them in a particular gauge. And given that multiple effective met¬ 
rics can satisfy the same governing equations, one must realize that under 
a gauge transformation, one is referring to the transformation of one’s par¬ 
ticular choice of effective metric. Furthermore, there is the potentially more 
dangerous caveat that the entire analysis of this section assumes smooth 
transformations. Under a transformation that is singular on the worldline. 


Rather than saying the equation of motion (351 is invariant, self-force literature usually 


talks about a transformation of the self-force. At first order, the equation in the new gauge 
is -t Ray/sUo The force is given by 


= Ff - 


d-To 


i± 


I TfU, a /-u p 

H“ -ti ai/'/3^0 Cl-L^O 


where the second equality follows from With z 


'M _ 


(138) 




M 

1 ± 


the left-hand side of the equation of motion, the equation’s invariance is transparent. 
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the singularity structure of the field can be significantly altered, the trans¬ 


formation laws for the deviation vectors are no longer given by Eq. (133), 


and the relation between the “nice” effective metrics in the two gauges will 


not be given by Eq. (135). 


6.3 Gauge in the self-consistent approximation 

I now consider the freedom in the self-consistent expansion, which is slightly 
thornier than the Gralla-Wald case. 


6.3.1 Transformation of the metric and the worldline 

Eirst, I note that the self-consistently determined worldline transforms 


according to Eq. (132). Unlike in the Gralla-Wald case, neither nor z'^ 
are expanded around Zq. 

Now to the metric. Let us adopt the passive view. Suppose we performed 
the self-consistent expansion in two slightly different coordinate systems 
and x'^. In the unprimed coordinates, the expansion is performed by writing 
as gfj.u{x, e; z) and expanding for small e while holding x^ and 2 ^ fixed. 
In the primed coordinates, we write g'^^ as g^^(x',e;z') and expand while 
holding x'^ and z'^ fixed. In terms of the perturbations = g^y — and 
= &'iiu ~ have expansions 


hi,y{x, e; z) = ^ z), 

h'^y{x, e; z') = ^ z') 


(140) 

(141) 


The metric in the two coordinate systems are related, as in Eq. ( 128| ), ac¬ 
cording to 

f^rpOL 

= ^^go.M^^e;ziz% (142a) 

= g^luix', e; z{z')) + C^g^y{x\ e; z{z')) 

+ 2^\9^^u{x' 1 e; z{z')) + O(e^). (142b) 

Expanding gfj,y{x', e; z{z')) returns 

^fiu^x , c) c \p'fii/(,x , z') ^y{x )] 4“ c 5 ^) 4~ ^^29^u{x ) 

\^\i9^iv{x')\ -l-O(e^). (143) 


In Eq. (143) I have simply left z^{z') as z^. Just as in the Gralla-Wald 
case, the term C^^hj^y{x'; z) introduces a mass dipole moment relative to 
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the reference worldline—in this case, relative to z^. But if we now expand 
z^{z') around the perturbation reads 


h'^y{x\t) = e[h]^^{x'\z') ++ e^[hl^{x'■, z') + C^^g^^{x') 

+ (Ai + '^ii)^]iv{x']z') + + O(e^). (144) 


where ] z') = 5hj^i,{x'; z 'uses the Lie derivative introduced in 

Sec.j^ The two Lie derivatives + £^^)h^^^{x'] z') do not precisely cancel 
one another, but £^^ does precisely remove the mass dipole moment induced 
by simply put, moves everything relative to the worldline, and £^^ 
moves the worldline by the same amount, such that there is no net change 
relative to the worldline. In other words, as we expect, contains no 
mass dipole term when written as a functional of the transformed worldline 
z'^. The fact that the dipole moment cancels in this way can be seen more 
explicitly by acting with + £^^) on the leading term in ~ 


2m5^ 
la :*—^*1 ’ 


it can also be inferred from the result (216). 


Equation (144) is an expansion at fixed x'^ and z'^. Hence, we can 
read off the coefficients of e"" and say that the individual terms in the two 
expansions (140) and (141) are related as 


h'luix'; z') = hp^{x'] z) + C^^g^y{x'), (145) 

z') = hl^{x'; z) + (£^, + £i^^)h\^{x'-, z') 

+ ^^29^iy{x') + ]^Cl^g^u{x'). (146) 


We can see from this analysis that a gauge transformation acts quite 
differently here than in the Gralla-Wald case. In the Gralla-Wald approx¬ 
imation, the zeroth-order worldline, on which the singular field diverges, 
is invariant; in the self-consistent approximation, a gauge transformation 
actually shifts the curve on which the singular field diverges. 


6.3.2 Transformation of the singular and regular fields 

As in the Gralla-Wald case, I define the split into singular and regular fields 
in the primed gauge in the simplest way that preserves the properties of 
the split. In other words, we must have the effective metric g^i^^ transform 
as any ordinary smooth tensor field would: g'^^^x') 
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which implies g'f,^{x', e; z') = ^^gaj 3 {x{x'),e] z{z')), and from there, 


z') = h^l{x'\z') + C^^gf,y{x'), (147) 

z') = h^l{x'] z') + (£ 5 , + £^^)h^l{x'-, z) 


At the same time we must satisfy Eqs. (145)-(146), which leaves the singular 
field to transform as 


h%{x'-, z') = /), (149) 

= hll{x';z') + (£g, + £^,)hll{x'-,z'). (150) 

With these choices, the worldline z'^ is a geodesic of and g'^^ is a 
smooth solution to the vacuum Einstein equation. 


6.3.3 Governing equations in alternative gauges 


Again as in the Gralla-Wald case, by design, the transformation laws (147)- 


(150) ensure that the governing equations of the self-consistent approxima¬ 


tion are invariant under a gauge transformation: in any gauge A", the regular 
field satisfies 6R^^[h^^\ = 0 , 6R^^,[h^‘^] = where is 

the linearized Ricci tensor in the X gauge; the center-of-mass worldline on 


which the field diverges satisfies (28); and assuming well-posedness, the full 


field off the worldline can be computed using the puncture scheme (98)- 


(100), with the replacement —)• —26R^^. The only change is to the form 

of the second-order puncture. The nature of this change can be inferred 


from the analogue of Eq. (139b), which reads 


6Rf,4{C^, + = 87t{C^, + %)r;, - 25^R^,[h^\Ah^^]. (151) 


SI 


,Rli 


This is just as in the Gralla-Wald case but for the presence of the £ term. 
The two Lie derivatives induce a change in the skeletal stress-energy, and 
hence in the seed solutions contained in the puncture, but unlike in the 
Gralla-Wald case, this change does not alter the mass dipole moment; the 
presence of ensures the mass dipole moment remains zero, as discussed 
above. So /i^®®'^(x; z, M) vanishes in all smoothly related gauges, and only 
h^if^'^{x', z,6m) is altered. 
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6.4 Gauge in the osculating-geodesics approximation 


Gauge freedom in the osculating-geodesics approximation is more compli¬ 
cated than in either of the previous two cases. Here, there is freedom at two 
levels: in the calculations in the Gralla-Wald expansion at each osculation 
instant, and in the calculation of the self-consistent perturbation h^y{x]z) 
that stitches together the Gralla-Wald expansions. 

If we always use the same gauge at the two levels, then there is no 
real complication: a gauge transformation at the self-consistent level can be 
expanded out to find the induced gauge transformation at the Gralla-Wald 
level in a straightforward way. However, in principle, calculations at the 
two levels can be performed in different gauges. In that case, one must have 
good control over the relationships between the gauges, since the final self- 
consistent approximation z) = Yln=i satisfies the Einstein 

equation through order only if z^ is in the same gauge as 

Suppose we wish to calculate h'^^{x‘, z) in a gauge X (denoted by un¬ 
primed symbols), but at an osculation instant r we wish to calculate the 
Gralla-Wald perturbations Zo(r)) in n gauge Y (denoted by primed 

symbols). To find the curve z^ in gauge X, we need to solve 


D^z^^ 


= z ) = eFf (r; Zo(r)) + ("t; 2:o(r)) 


(152) 


with the forces in gauge X. But as input we have only the T-gauge per¬ 
turbations h'J^^{x; 2 o(t))- For simplicity, consider only the first-order force in 
gauge X. It is related to T-gauge quantities as [see Eq. (138)] 


Fiij) = F[^{t) + 


(153) 


where = -\Po''[2h'^a-AF ZQ(r)) - To know the 

correct force for our self-consistent evolution, we need to know the gauge 
vector ^2*. However, this gauge vector is constrained by the osculation con¬ 
dition: we must have z^{t) = and n^(r) = ttp^^^(r). Hence 

should impose (z(r)) = 0, reducing Eq. (153) to 


we 


Ff = Ff -F 


dro 


(154) 


Using this relation, one can calculate the force that appears in Eq. (152) 
using the force calculated from the Gralla-Wald expansion in gauge Y —but 


one 


must know at the osculation instant r. The same analysis can be 


straightforwardly generalized to the second-order force F^(r; 2 ;o(t-))- 
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7 Equations of motion from a rest gauge 

In the algorithm of Sec. the nth-order equation of motion is obtained by 
solving the + l)th-order Einstein equation in the buffer region. However, 
at least in certain cases, we can instead obtain an equation of motion using 
only the nth-order solution. This is accomplished by finding the relation¬ 
ship between one’s desired gauge—call it the practical gauge —and what I 
call a rest gauge, which is, essentially, any of the gauges that have been used 
to study tidally perturbed objects in the context of matched expansions 
(see Ref. and the many references therein for examples of these stud¬ 
ies). Methods along these lines were the first ones used to derive the MiSa- 
TaQuWa equation [23] and currently, they are the only ones that have been 
used to derive second-order equations of motion [161 US [211 EE] , although 
different notations and descriptions have obscured the shared underlying 
features of these derivations. 

Throughout this section, I specialize to the case of an approximately 
spherical, nonspinning object whose leading-order dipole and quadrupole 
moments vanish: = 0. I only sketch the calculations in this 

section refer the reader to Refs. [I8l|22| for the details of the calculations in 
this section. 


7.1 Metric in a rest gauge 


By a rest gauge, I mean a gauge in which the object is manifestly at rest 
relative to a worldline and manifestly centered on that worldline. To 


motivate this idea, I return to the themes of the introduction. In Sec. 1.1.2 


I argued that any equation of motion can be written as the geodesic equa¬ 
tion in some smooth piece of the metric. In Sec. |1.3[ I recalled the results 
of Thorne and Hartle, who showed something slightly stronger: that the 
equations of motion and precession for any compact object can be written 
as those for a test body immersed in some effectively external metric. In 


Eq. (82), the self-force program has recovered their result for the center-of- 
mass motion, and it has shown that the “external” metric in their formalism 
is in fact (at least through order er) the effective metric g^j^p = g^y + h^y. 

At second order in e and beyond, the algorithm of Sec. has not yet 
revealed how g^j^y relates to Thorne and Hartle’s external metric. But if 
we can establish the relationship, we can establish the equation of motion. 
Notably, they derived their equations by working in an inertial frame of their 
external metric; if the object possesses no spin or quadrupole moments, then 
there is no force on the object in this frame, and it is what I call the rest 
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gauge. Why I call it a gauge rather than a frame will become clear as we 
move forward. 

The metric in this gauge is most easily constructed through the inner 


expansion (39) rather than the outer. To do so, I first note that just as the 


matching condition fixes a lowest allowed power of r in the outer expansion, 
so it determines a highest power of r in the inner expansion: expanding 
the nth-order inner perturbation in the buffer region, we find 

e"'J2p{e/r)PHll^u{t,n^); the condition that no negative powers of e occur in 
the outer expansion then determines 




p<n 




(155) 


Here and in what follows, I assume some locally Cartesian coordinates cen¬ 
tered on the center-of-mass worldline z^. Given an appropriate choice of rest 
gauge, the coordinates will later be identified with Fermi-Walker coordinates 
centered on z^. 

For our approximately spherical, approximately nonspinning object, the 
inner background metric, when expanded in the buffer region according to 
Eq. (42), is equal to the Schwarzschild metric through order 1/f^. It has 


the schematic form 


g 


,obj 


m m? ,_a. 

ri+ — + ^ + ^ + r% 


(156) 


here I do not specify any particular coordinate system, though I insist that 
the coordinates are mass-centered, such that no mass dipole moment appears 
in the expansion of every term in Eq. (156), at the orders displayed, is 


fully determined by m. Using our previous results, we know m is constant. 
Hence, we can find the perturbations in the buffer region using the 
well-developed formalism of perturbations of Schwarzschild. Furthermore, 
given the assumption that the perturbations are quasistationary (because, 
as explained in Sec. the inner expansion scales space but not time), time 
derivatives appear as higher-order terms in the Einstein equations. That is, 
one first solves the time-independent linearized Einstein equation 


6R%[g°^KH^]=0 


(157) 


for with t fixed; here 6Rfj,u[g°^^, H^] = 0 means the Ricci tensor lin¬ 
earized off the background and the superscript “0” means “terms in 
6Rnu containing no t derivatives”. Next, one solves 




( 158 ) 
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for at fixed t and for the time derivative of where the superscript 
“1” in srU means “terms in (5-R^i/ containing one t derivative”. At this 
point, the pattern should be obvious. 

Just as when constructing the outer expansion in the buffer region, we 
here solve a sequence of spatial differential equations. Furthermore, it is 
well known that in an appropriate gauge, stationary vacuum perturbations 
of Schwarzschild behave as or at large f, where i is the spherical 

harmonic index, just as in the outer expansion. Let us label the 
solutions as and the solutions as Combining 

this with Eq. (155), we find that at first order, where there is no source, the 


solution behaves as 




+ 0(l/f3) 


(159) 


and 


It is well known that the monopole and dipole solutions 
are pure gauge. Hence, I set them to zero. The gauge-invariant content in 
is a correction to m, and in it is a perturbation toward 

the Kerr metric—that is, a small spin. With a bit of work, one can show 


from Eq. (158) that these perturbations are independent of t, and we can 


freely set them to zero, absorbing the mass correction into m and specifying 
that the object remains unspinning at this order. Physically, this time- 
independence corresponds to the fact that physical changes in mass and 
angular momentum are caused by tidal heating and torquing, which are 
caused by nonlinear effects at much higher order in e. So in the end, we 
have 


= 0 . 


(160) 


Since H}„, vanishes, Eq. (158) is again the linearized vacuum equation: 




= 0. The solution behaves as 


~ + fOH'2,o,+ ^ 1^2,0,- ^ 0(l/f2). (161) 

r 


Again, 


r2,0,- 


Lf^i, can be gauged away. The invariant part of 


2 , 0 ,- 


and 

is another correction to m, which again can be absorbed into a redefinition 
of m. Writing ■, we then have 






{t, r)h^^ -|- 0{\/r^). 


(162) 


lim 


f^oo {t,r) can be written in terms of two STE functions £ij{t) 

and and in the expansion of {t,f) for large r, every term is 
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directly proportional to one of those two functions (i.e., constructed from 
contractions of one of them with 6ij or €ijk)- These two quantities can be 
interpreted as quadrupole tidal moments in the neighbourhood of the object, 
and from the matching condition they will be identihed with those of the 


external background metric, defined in Eq. (46). 


The third-order Einstein equation now reads 

Eollowing the same steps as at the previous orders, we find 


(163) 


^f^H' ... 


+ f^H, 


2,2,H 


.’. +r (t, r)h*^ + C>(l/r). 

(164) 

r) can be written in terms of two STE functions £ijk{t) and Bijk{t), 
such that every term in its large-f expansion is directly proportional to one 
or the other. These two quantities can be interpreted as octupole tidal mo¬ 
ments in the neighbourhood of the object, and they will be identified with 
those of the external background metric, which are constructed from covari¬ 


ant derivatives of the Riemann tensor of g^j^y. Equation (164) additionally 
contains corrections to the quadrupole fields, both in the form of the time 
derivatives tij and Bij appearing in , and in the field which 




can be written in terms of two STE functions, call them 5£ij and 


like 

Now let us put the results together. After we rewrite it in terms of 
unsealed coordinates r = ef and re-expand in e, the metric g^y = g^^ (f) -|- 
e^-f^L(^) + + O(e^) reads g^,y = g^,y{r) -h eh'^y{r) + e^h'^^{r) -h 

e^h'^y{r) + O(e^), with 


(165) 


where primes refer to quantities in the rest gauge, and g = diag(—1,1,1,1). 
‘-b’ signs here mean “plus terms directly proportional to”. Eor compactness, 
I have omitted the magnetic-type moments Bij, Bij, 6Bij, and Bijk', they ap¬ 
pear in exact analogy with the electric type moments £ij, £ij , 5£ij, and Sijk ■ 
All terms in the tableau, at the displayed orders, are directly proportional 
to one of these moments and/or to a power of m. Every term in the hrst col¬ 
umn is £ = 0, every term in the second is ^ = 2, and every term in the third 
is either i = 2 (the £ij and Bij terms) or .^ = 3 (the Sijk and Bijk terms). In 


9 = 

g 

-b 

r^Sij 

+ 

^ i^ijk “1“ ^ij) 

+0{r^ 


m 

-b 

mrSij + r‘^6£ij 

+ 


+0{r^] 

to 

m? 

7T 

-b 

m?£ij + mr6Sij 

-b 

rn^r{£ijk -\-£ij) 

+0{r^] 

h'^ ~ 

rrS 

y3" 

-b 

^Sij + mHS^, 

-b 

'kkl {^£ijk “1“ 

+0{r) 
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this gauge, the object is manifestly at rest on z^\ there are no mass dipole 
moment terms, which have the even-parity dipolar behavior n'‘ jr^^ nor any 
acceleration terms, which have the even-parity dipolar behavior rre*. 

By choosing the gauge of the inner expansion appropriately, we can put 
the external background metric given by the top row of Eq. (165), in 


precisely the form (45) that it takes in Fermi-Walker coordinates—except 


that no acceleration terms appear. Since we know that the center-of-mass 
worldline is accelerated in this tells us that Eq. (165) is actually in¬ 
complete. The construction we have followed here is valid for any worldline 
, regardless of whether it is accelerating in g^^ or not. The fact that 
we see no explicit acceleration terms is equivalent to the empty statement 
that any worldline z^ can be written as a geodesic in some smooth piece 
of the full metric. Ergo, at this stage we know nothing at all about z^, 
and to make g^u the external background spacetime, we must insert into it 
by hand the acceleration terms appearing in Eq. (45). At the same time. 


in our choice of gauge, no acceleration terms appear in the full metric 
Therefore, we must insert Fn terms into the perturbations such that the 
sum g^v{r) -)- eh'^^{r) + c^h'^y{r) + e^h'^^{r) + 0{e^) contains no acceleration 


terms when one applies the expansion 

The form of g^u in the rest gauge invites us to make a new split into an 


effectively external metric = g^^ + that contains the m-independent 


fLU 

tidal terms in Eq. (165), plus a self-field = gf^y — Q^n,, every term of which 


is directly proportional to a power of m. With an appropriate choice of rest 
gauge, looks exactly like Eq. (45) with none of the acceleration terms; 


that is, our gauge makes the Fermi coordinates refer to proper distances, 
times, Riemann curvature, and parallel transport defined with respect to 
Qfiv We have 


Sit = -1 - + C>(r^), 

Qta = O(r^), 

gab = 5ah - l^aibjX'x^ + 0{r^), 


(166a) 

(166b) 

(166c) 


where = Ri^aup + edRfj_aui3[^^^] + C>(e^) is the Riemann tensor of 

As per the discussion of acceleration terms in the preceding paragraph, the 
new “regular field” reads 


= 2Fini - 6Roioj[i)^]Fx^ + ©(r^), 
flto = -l5Roiaj[i)^]x^x^ + 0{r^), 

= -ISRaibjlh^Wx^ + 0{r^), 


(167a) 

(167b) 

(167c) 
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+ Oir^), = 0{r^), 


(168) 


and = 0{r). For visual clarity, I have truncated these equations at 


lower orders in r than Eq. (165). 


In summary, in a rest gauge the metric perturbations take the form 






Rn 


(169) 


where is made up entirely of terms containing explicit factors of m”, 
no mass monopole corrections, no spin moment, no quadrupole moment, 
no mass dipole moment, and no acceleration terms, and is given by 


Eq. (167)“(168). Equation (169) is the general solution in the buffer region 


outside an approximately spherical and nonspinning object, just written in a 
particular gauge. In this gauge, is manifestly a geodesic of an effectively 
external metric, and is being defined as that effectively external metric 
in which is a geodesicj^ A priori, need not be related to the 
latter was defined in a very different manner and is not yet known to be the 
metric in which the motion is geodesic (except at zeroth and first order). 


7.2 Self-consistent approximation: worldline-preserving gauge 
transformations 

Metrics in local rest gauges are common in the literature, usually derived in 
the context of tidally perturbed compact objects. The goal of the self-force 
game is not simply to construct such a metric, but to find the equation of 
motion in whichever “practical gauge” one wants to use to compute a global 
solution. 

Let us take the Lorenz gauge as the practical gauge. I wish to find 


a unique gauge transformation between the rest-gauge perturbations (169) 


and the Lorenz-gauge perturbations (70)-(73) 
and ^2 satisfying 


That is, I seek gauge vectors 








1 

pv • 


(170) 

(171) 


^®Note that without some additional input beyond that definition, the split of into 
g^^u, and is quite ambiguous. Suppose there were a zeroth-order force acting on 
the object. One could still write the equation of motion as a geodesic in some smooth 
piece of the metric, but to do so, one would have to shift part of into one would 
not simply be splitting the perturbations h'fi,. With the present setup, the ambiguity 


that Fq = 0. 


is lifted by assuming the expansion (251 and utilizing the independently determined fact 
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Note that here I do not utilize the third-order perturbation, though its form 
was found in the rest-frame gauge. What is its relevance? Only that 
it contains no mass dipole moment term; in other words, it shows that 
with this choice of gauge, the correction to the mass dipole moment can 
unambiguously be set to zero to define a center-of-mass-worldline. 

These equations are to be solved subject to one crucial condition: they 
must preserve the location of the worldline. The rest-gauge solution shows 
that the object is at rest on some worldline, the Lorenz-gauge solution is 
written in coordinates centered on some worldline, and in order to claim 
that the worldline of the Lorenz-gauge solution is the desired center-of-mass 
worldline, we must enforce that it be identical to the worldline of the rest- 
gauge solution. In the simplest scenario, Eqs. (170)-(171) can be solved 
with a smooth gauge transformation of the form 




(172) 


p>0 £<p 


Given Eq. (132), the condition that the worldline not be altered is that 


Cooit)=^± = 0. 


(173) 


I call a transformation satisfying this condition a worldline-preserving trans¬ 
formation. This condition was first used in Ref. [23]. A more general con¬ 
dition is used in Refs. dSlIIE]. 


By working through sequential orders of r in Eqs. (170)-(171), one 
uniquely determines ^((^^(t) (up to residual gauge freedom in the rest gauge 
and Lorenz gauge). I do not present details of that process here, but the 
explicit solution to the first-order equation (170), with particular choices 
of rest gauge, can be found in various references; see, e.g., Ref. [Hj. The 
explicit solution at second order is presented in Ref. [T8|. What is the es¬ 
sential result in these solutions? At first order, they uniquely determine 


Ff, which appears via Eq. (167), to be the MiSaTaQuWa force (66) (with 
S'* = 0), as we already know. At second-order, they uniquely determine the 


second-order force Ft, which appears via Eq. (168), to be 


Ff = \dihu‘^ - dth^‘^ - 


(174) 


Eollowing Appendix A.l, the results for F(* and can be combined to 
write the equation of motion = eF^ e^Flf + 0{€^) as the geodesic 

equation = 0{e^) in = g^y + h^y. This is the result promised way 
back in Sec. [TXTl 
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What is the meaning of this result? In the rest gauge, we knew that 
was a geodesic in Now we know it is a geodesic in Either of 
or can be thought of as a “nice” effectively external metric. Are they 
the same metric? Not necessarily. The two geodesic equations only restrict 
the metrics and g^iy on the worldline and their first derivatives on the 
worldline; they place no constraint on higher derivatives. But one of the 
results of the process of finding and ^2 is that the two metrics are related 
by 


+ 0{r‘^) (175) 

+ C>{r^). (176) 


In other words, they are equivalent up to possible O(r^) differences. Are 
they the same at order and higher? The answer depends on precisely 
how the rest gauge, and the tidal moments 5£ij, 5Bij (and at higher order 
in r, octupolar and higher order moments) appearing in is constructed; 
this construction is, unfortunately, not unique. Reference [18] demonstrates 
explicitly that the two effective fields can differ at order r^. 

Is there a way to determine which effective field is the “best one”? As 
described in Sec. 3.5, the most promising route is by solving the field equa¬ 


tions outside an object with spin and higher moments, since those moments 
will couple to some tidal moments, and we can determine which metric the 
tidal moments belong to, and thence the metric in which the object moves 
as a test body. But just with the results at hand, we can confidently say 
that g^y has all the nice properties one might desire of it for a spherical, 
nonspinning object: it is a vacuum solution, causal on the worldline, and 
the object moves on a geodesic of its geometry. 


7.3 Gralla-Wald approximation 

Starting from the self-consistent results, one can readily derive the second- 


order equation of motion (36) following the expansion procedure of Ap¬ 
pendix |A.2[ 

Alternatively, one can derive a second-order equation of motion directly, 
as was done by Gralla [22j. In the Gralla-Wald case, the rest gauge is 
constructed in coordinates centered not on but on Zq. The results are 
the same as in Sec. 7.1, but all acceleration terms are set to zero in g^^y (and 
likewise for the Fn terms in that cancel them). When transforming to 
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a “practical gauge”, one uses 


(177) 

(178) 






and one does not impose the worldline-preserving condition. The gauge 
transformation is allowed to move the object relative to Zq. 


Z 2 = 0 in the rest gauge, we have from Eq. (|l33l) that in the practical gauge 


Since z'j^ = 


(s) = -^r(^o), (179) 

4(s) = -^2^(^o) + \ei{z,)d,eA^o), (180) 


and the covariant second deviation is Z 2 p{s) = —^2 (^ 0 ) + i/ii{zo). 

Following Gralla, one can use this procedure to obtain results in any 
gauge smoothly related to a particular rest gauge. Defining regular fields 


= (181) 
+ ^i29tMu + (182) 

one finds evolution equations for Zn in terms of by taking appropriate 
derivatives of these relations; any desired gauge condition is imposed on 
h^y is computed via a puncture scheme, and Zn is computed from h^y. 
Because of the particular method of construction, in which the rest gauge 
is defined with respect to Zq rather than z^, the evolution equations for Zn 
do not correspond to an expansion of the geodesic equation in the metric 
Q/iiy = 9/112 + ll^^-l Gralla defines; this lack of geodesic motion in g^y + \}^y 
should amount to a slightly different definition of \)^y than the one used in 
the self-consistent expansion, although no one has yet carried out a detailed 
comparison of Gralla’s second-order results to those of the self-consistent 
method. 


8 Conclusion 

8.1 Summary 

In this review, we began with the idea of a point particle interacting with 
its own backscattered field; we ended with the idea of using laws of gauge 
transformation to obtain laws of motion for extended objects. And yet, 
reassuringly, from the latter description we have recovered the former. 
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Along the circuitous way, I have emphasized several core ideas. First, the 
laws of motion of an extended object, to all orders in perturbation theory, 
are determinable from the metric outside the object; all information about 
the object’s shape and internal composition is encoded in a discrete set of 
multipole moments, themselves defined from the form of the metric in the 
buffer region outside the object. To obtain an equation of motion of order 
in perturbation theory, moments of multipole order i = n are required. 
Second, the laws of motion found in this way show that at least through sec¬ 
ond order in perturbation theory, a small, sufficiently spherical, sufficiently 
slowly spinning compact object moves on a geodesic of an effectively exter¬ 
nal metric that (i) satisfies the vacuum Einstein equations and (ii) behaves 
causally on the object’s representative worldline. Furthermore, the analysis 
of the metric outside the object provides a clean way of separating the met¬ 
ric into two constituents: a suitable effectively external metric satisfying the 
above properties, and a self-field that, loosely speaking, locally characterizes 
the object. 

This treatment weds several traditions that run through the history of 
the problem of motion in general relativity, from the characterization of mo¬ 
tion of a material body in terms of suitably defined multipole moments in the 
tradition of Mathisson |37] and Dixon |36j , to the method of algorithmically 
characterizing the metric in terms of multipole moments in the multipolar 
post-Minkowskian theory of Blanchet and Damour ISIEH], to the derivation 
of laws of motion of asymptotically small objects from the Einstein equation 
outside the objects in the tradition of Einstein [21], D’Eath |32|) and Thorne 
and Hartle |l8| . The perturbative treatment here also complements the non- 
perturbative treatment of Harte [38]. Just as Harte brought Dixon’s project 
to fruition by Ending a utile separation of the metric into a self-field and an 
effectively external metric, so the same can be said about the method here 
bringing to fruition the project of Thorne and Hartle. 

An important aspect of any perturbative treatment of motion is the 
relationship between that motion and gauge freedom. It is well known that 
on short timescales, the self-force is pure gauge; it can be freely set to 
zero by a suitable choice of gauge. Equivalently, one can say that on short 
timescales, the perturbed worldline of the object can be transformed to a 
geodesic of the background spacetime. 

What physical role, then, does the self-force play? On short timescales, 
one can say the following: the deviation from geodesic motion appears ex¬ 
plicitly in the second-order field. Hence, in order to get gauge-invariant 
information about the second-order field from the solution to the Einstein 
equation in any given gauge, one must include the effect of the self-force. 
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(Note that the argument doesn’t run in reverse. Since the first-order devia¬ 
tion is determined entirely by the hrst-order gauge choice, one does not need 
the entire second order perturbation to obtain invariant information involv¬ 
ing the first-order deviation m-) However, on these short timescales, the 
perturbative correction to the motion need not be accounted for to obtain 
physical information from the first-order metric. 

As I have argued, this situation changes when one considers long scales. 
If one seeks a solution valid on a large domain, such as over the course of a 
binary inspiral, then one must account for the correction to the motion in 
order to obtain a well-dehned perturbation theory. If one were to attempt 
to put the perturbative shift in the motion into the second-order field, then 
on this timescale the second-order held would grow larger than the hrst- 
order held. Hence, one must incorporate the effect on the motion into the 
hrst-order perturbation. In this sense, the perturbed worldline itself is quasi¬ 
invariant under a gauge transformation on long timescales; it may only be 
shifted by a small amount of order e, while the shift due to the self-force is 
a much larger effect, of order ~ I/e (or ~ e^, depending on whether one is 
looking at orbital phase or orbital radius in an inspiral, for example). 

For this same reason, if one is interested in long timescales, one need not 
include the second-order perturbation in order to obtain gauge-invariant 
information about an inspiral. For any given short patch of that inspiral, 
the effects of the self-force will be pure gauge; but the accumulation of those 
effects over long timescales is invariant within the class of gauges that are 
well behaved on these timescales. 

8.2 Future directions 

Although the foundations seem in place and the second-order results seem 
sufficient for the practical purposes we can presently imagine, several open 
avenues remain to be explored. 

First, at present we have obtained second-order equations only for ob¬ 
jects whose spin and quadrupole moments vanish at leading order. It would 
be worthwhile, and astrophysically relevant, to obtain second-order equa¬ 
tions of motion for more generic compact objects. As is well known [36], 
the object’s leading order quadrupole moment would generate a force by 
coupling to the external curvature. Additionally, there may be an effect due 
the correction (55* to the object’s spin. 

The mention of this subleading spin raises another question to consider: 
what is the physical content of the perturbed multipole moments? I have 
discussed above how the monopole correction is pure gauge (at least on 
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short timescales). Will higher corrections, or corrections to higher moments, 
contain more physical information than this? We can expect that they will: 
tidal heating and torquing, for example, should create physical corrections 
to the mass and spin |80l l8T] , and at least for a material body, the external 
tidal fields should correct the higher moments by physically deforming the 
object [HSHHlj. Ultimately, to model the motion of a particular class of 
objects, such as realistic neutron stars, one will need to explicitly match 
an inner expansion to the outer expansion I have discussed; the matching 
procedure will uniquely identify the multipole moments of the objects as 
they appear in the metric in the buffer region. To this end, one might 
consider metrics of tidally perturbed neutron stars, such as those described 
in Refs. [82l[83]. 

One might also try to describe less compact astrophysical objects, or 
more exotic objects, by broadening the scope of the perturbative expansion. 
Rather than assuming the object is compact, such that its linear dimension 
d is of the same order as its mass m, one might consider a two-parameter 
family corresponding to m and d. This would alter the orders at which 
various multipole moments appear in the buffer region, since they scale as 
md^. One could even try to describe approximately string-like objects by 
examining a limit in which two of the object’s linear dimensions go to zero 
while one linear dimension remains finite. 

To gain further insight into the physical content of the multipole mo¬ 
ments, one could also relate them to the moments appearing in other for¬ 
malisms, such as the non-perturbative formalism of Harte j38| . Another 
point of comparison with Harte would be the definition of the effectively 
external (i.e., the “regular”) field. As I have discussed, the choice of ef¬ 
fectively external field is far from unique, and in fact an infinite number 
of choices could be made that would still guarantee the center-of-mass mo¬ 
tion is geodesic in the effective metric. Here, I have taken as my guideline 
that the effective metric should be a vacuum solution that is causal on the 
worldline and in which the motion is geodesic; these conditions still do not 
uniquely identify an effective field, but they narrow the range of options, 
and I have shown how a choice satisfying these conditions arises naturally 
in the process of solving the Einstein equation in the buffer region outside 
the object. However, Harte makes a different choice, which does not satisfy 
the vacuum Einstein equation. How are the two related? The answer is not 
obvious. One way of establishing a stronger relationship would be to find 
perturbative equations of motion for objects with higher multipole moments: 
because moments of higher i couple to higher derivatives of the metric, they 
feel much more of the field than does a monopole. With Harte’s definition 
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of effective field, the equation of motion is precisely that of a test body in 
the effective metric. Would the same be true in the effective metric I have 
defined here? 

In addition to comparison with Harte, the results I have described should 
be related to Gralla’s more closely related, perturbative results [22]. At 
second order, Gralla makes a choice of regular field that is a vacuum solution 
(at least to the order in r to which he defines it), but in which the center- 
of-mass motion is not geodesic. 

Although I have emphasized the idea of splitting the physical metric into 
a self-field and an “effectively external” metric, and although I have linked 
that split to a generalized equivalence principle, the array of choices of reg¬ 
ular field should make clear that there is a danger of over-interpreting the 
physical meaning of any particular choice, no matter how nice its proper¬ 
ties. This risk is also present because in general, a regular field satisfying 
nice properties on the object’s worldline is acausal when evaluated away 
from the worldline. As a practical matter, the freedom to define different 
effective metrics may make future comparisons of self-force results to post- 
Newtonian results more hairy; given the vast freedom—even while main¬ 
taining the property that the motion is geodesic in the effective metric, for 
example—it is remarkable and fortuitous that agreement has so far been 
found for so many different effects that appear to rely on particular choices 
of this field jSSHSZj- 
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A Expansions of the geodesic equation 

A.l Expansion in powers of a metric pertnrbation 

In this appendix, I examine the expansion of the geodesic equation in any 
sufficiently smooth metric g^i^; the treatment is generic, not specialized to a 
spacetime containing a small object. I expand only in powers of a sufficiently 
smooth metric perturbation; I do not expand the worldline itself. Hence, the 
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analysis is meant to apply to the self-consistent representation of motion, 
not to the Gralla-Wald representation. 

The geodesic equation reads 


(j 


(183) 


where s is a potentially non-affine parameter on the curve z^{s), 

is its tangent vector field, is the Christoffel symbol corresponding to 

gpu, and K = £ In 

If we now write the metric as the sum of two pieces, = g^i, + hpy, 
and if we take s = t, the proper time on z^ in g^iy, and if we rewrite the 
geodesic equation in terms of covariant derivatives compatible with g^^, we 
find 

= -AT^u'^uP + (184) 

where and 


Ar?, = - a = (2/.W;,) - 


(185) 


is the difference between the Christoffel symbol associated with the full 
metric gpp and that associated with the background metric g^i,. With r as 
a parameter, n becomes 


Y^l - hpuUf^u’^ 


(186) 


So far no approximation has been made; Eq. (184) is exact. If we now 
expand AT^p and k in powers of hp^, we find 


a 


a 


= ~ ^"'^)( 2 ^ 5 (/ 3 ; 7 ) - h0^.s)u^v? - ^h^^-su°‘u^u^u^ 

— ^hpiyhi3.y-su°‘u^u^u^u^u‘' — h^^u°‘a^v? + 0{h^). 


(187) 


This equation is complicated by the fact that the acceleration appears on 
both sides in a nontrivial way. To disentangle the acceleration from the 
perturbation, I assume that a^, too, has an expansion in powers of hp^, 

«'' = <n + «quad + ^(/'^)> (1^8) 
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where is linear in and aqu^d quadratic in it. Substituting this 
expansion into Eq. (187), one finds 


(189) 

(190) 


«Un = ~ 

“quad = “ V;5)u^u'^, 

where Summing these, we have 

- h\){2hs[i3-^) - hi3^.s)u^u'^ + 0{h^). (191) 

As applied to the case of the effective metric = g^y + (i.e., 

replacing /i^i/ with /ij)'^), Eq. (191) agrees with the second-order self-forced 
equation of motion derived in the body of the paper. 


A.2 Expansion in powers of a metric perturbation and a 
worldline deviation 

In the last section, I expanded the geodesic equation while holding the so¬ 
lution Ze (s) to that equation fixed. I now expand Ze (s) as well. This proce¬ 
dure yields a sequence of equations for the terms in the expansion of Ze{s), 
suitable for a Gralla-Wald approximation. My approach to the expansion 
closely follows the treatment of geodesic deviation in Sec. 1.10 in Ref. j88| 

I first describe the geometry of the situation. Consider a family of world¬ 
lines 2 ;^(r, e), with each member Zeij) = z^{T,e) governed by the equation 
of motion (|184). Each member satisfies 




= F^{T,e), 


(192) 


where r is proper time on Ze, and is given by the right-hand side of 


Eq. (184). The family generates a two dimensional surface S with a tangent 


bundle spanned by and An important relation between 

these vector fields can be found from which implies CuV^ = 

0 = CyU^, and from there, 

vf^.yu’^ = u^^.yv^. (193) 

Now, we seek to describe the deviation of an accelerated worldline z^ir) 
from the zeroth order, geodesic worldline Zq(t) = z^(t,0). The first step is 
to expand the worldline in the power series 


2 ;^(r,e) = z^(t) + ez^ir) + e^z^ir) + O(e^), 


(194) 
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where 


(195) 


Zni^) = 

n! 


1 „ d'^z^ 


de^' 


(r, 0). 


We may also write the expansion as z^{T,e) = ^ :^e"'C^z^\zQ(r)- Note that 
here z^ is a scalar field equal to the /xth coordinate field on the surface S. 
The leading-order term is the family member Zq{t) = z^(r, 0). The second 
term is z'^{t) = CyZ^\zQ = v^{zo{t)), a vector on Zq. But at second or¬ 
der and beyond, a subtlety arises: unlike the first derivative along a curve, 
which is a tangent vector, second and higher derivatives are not immediately 
vectorial quantities. The function z^{T,e) describes a curve in a particular 
set of coordinates, and the corrections Zn depend on the coordinate system 
in which one defines z^{T,e). Since my notion of an object’s center is es¬ 
tablished with reference to a comoving normal coordinate system, I wish 
my covariant measure of the second-order deviation to agree, component by 
component, with 5 0) when evaluated in a normal coordinate system 


centered on Zq ; Sec. [^describes the utility of this choice when re-expanding a 
self-consistent approximation into Gralla-Wald or osculating-geodesics form. 
With that in mind, I define the vector 


1 Dv^ 


^ 2 


de 




and I seek an evolution equation for its restriction to Zq, 

4f(^) = W^“Uo(r)- 


(196) 


(197) 


Z 2 P is the second-order term in the expansion (194) when that expansion is 
performed in Fermi normal coordinates centered on Zq. 


In addition to the choice of coordinates, the expansion (194) depends on 


the particular choice of parametrization (r, e) of the surface S. A change of 
parametrization alters the direction of expansion away from Zq. Here, the 
parametrization is chosen such that r is proper time along each curve Ze(T), 
and a flow line generated by links points on different curves Ze (t) at the 
same value of r. When restricted to Zq, the parameter r is tq, the proper 
time on Zq. 

With all preliminaries established, I now proceed to find the evolution 
equations for Zq , z^, and z^. The leading term clearly satisfies 


dro 


= Ff^{T,0) = 0 . 


(198) 


For the others, I first find evolution equations for and and then 
evaluate the results on Zq. At first order, using Eqs. (192) and (193), we 
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have 


D‘^1 


dr^ 


= V 


;/3^ 




;7 


;7 




(199) 

( 200 ) 
( 201 ) 


where the second line follows from Eq. (193) and the third line follows from 
the Ricci identity and Eq. (192). Evaluating on Zq , I write this as 




= F^iro) - 


( 202 ) 


where E" = |.yp. This is a generalization from the usual geodesic devia¬ 

tion equation to the deviation between neighbouring accelerating worldlines; 
it is valid even if E^(r, 0) / 0. 

At second order, repeated use of Eq. (193) and Ricci’s identity leads to 
D^w°‘ 






(203) 




;7 


+ -R!^ ^pv-'y {v^u^v^v? — u^v^v^'v'^'^ 

- . (204) 


Evaluating on Zq and using Eq. (192), we can write this as 

= A“(to) - R^^y^pu [u^z^^Uq + 2u(‘2;fuo^ -h IR"^ y,py.^z{^u^’z'y 


(/^,,/3).,b„,7] 
«0 

(205) 


where F 2 = I 70 ^ Equation (205) describes the sec¬ 

ond deviation between neighbouring accelerating worldlines. In the case of 
neighbouring geodesics, it agrees with “Bazanski’s equation” in the form 
given in Eq. (5.9) of Ref. [89] . 

The quantities Fn appearing in Eqs. (202) and (205) can be straight¬ 
forwardly evaluated by performing the expansion h^y{x,e) = ehjj^^{x) + 
e^hj^j^{x) F 0{e^) in Eq. (191) and then taking covariant derivatives with 
respect to v^. The results are 


1 




(206) 




















and 


^2 — ~ 2^0 ~ UqUq — -Pq UqUq zf 

- (2^L;A - hlx-,u) 

+ P^h^ {2hl,.,x - hlx;p) «o«o- (207) 


As applied to the case of the effective metric g^u = g^iu + (i.e., 

replacing with Eqs. (202) and (205), together with Eqs. ( |206 ) 

and (207), are the second-order expansion of the motion that apply in a 
Gralla-Wald approximation. 


B Expansion of point-particle fields in powers of 
a worldline deviation 

In this appendix, I derive the linear terms in expansions of the point particle 
stress-energy Tf^(x;z) and the Lorenz-gauge retarded field /i^^(x; 2 :) when 
the worldline is expanded as z^{s, e) = Zq(s) + ez^(s) + O(e^), where s is an 
arbitrary parameter. I also establish the identification of these linear terms 
with the mass dipole moment terms found from the local analysis in Sec. 


B.l Stress-energy 

I write the stress-energy in the parametrization-invariant form [9] 



(208) 


where z) is a parallel propagator from the source point x' = z(s, e) to 


the held point x, and z^ = ^ 


Substituting the expansion (30) into this stress-energy tensor, we obtain 


T^^(x;z) = m f g^,{x, zo)g'^ 0 ,ix, zo)z^'z^' + ez'^'vf,>{g 2 ^g^^,z°‘'z^')\^=o 


'70 


X 6{x,zq) + eziVyi5{x,z)\e=o 

+ 0 {^). 


S' 


1 - e 


ZQS'zj Vyig 

Zq ZQf,' 


ds 


-ZqZOp' 


(209) 


In each instance, the evaluation at e = 0 occurs after taking the derivative. 
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I simplify this expression using the distributional identities V^/(5(x, z) = 
—g^,Vfj_5{x,z) and g^,.^,6{x, z) = 0 [H]. I also use the identity 




ds 


1 - za' 

^1 


( 210 ) 


which follows in the same manner as Eq. (193) 
The result of those simplifications is 


er"^(x; z) = eTf'"(x; zq) + e^5T^'^{x] zq, zi) + 0{e'^), 




a/3, 


( 211 ) 


with 


rf^(x; zo) = m I g^,g'^,u^'u^ 5{x, zo)dro, 


( 212 ) 


'70 

(5rf^(x; Zo, zi) = m f g2>g^, ^ - UqU^ uoynj ^ 6{x, zq) 

- 9Z'^'r^ix,zo)\dTo, 


(213) 


where u^{to) = and the parallel propagators are evaluated at (x, zo{tq)). 
I have simplified these expressions by reparametrizing Zq in terms of tq, the 
proper time on 70 , but note that this does not correspond to choosing the 
original parameter s = tq. Equations (211)-(213) are valid for any choice 


of parameter s, and Zi{to) actually depends on the original choice of s: a 
change of parametrization s —)• s'(s,e) will change the direction of z^, in 
particular changing whether or not z^ is orthogonal to Uq. 

To eliminate this dependence on the initial choice of parametrization, I 
rewrite in terms of the orthogonal part of z^, z^^ = {5^ + UqUou)zi. 

The result is 


= m [ 2u[)“ ufl6{x,zo) - UqU^ zJj_g^,Vy6{x, zq) dr^ (214) 


'70 


where u^j_ = ^he part of z^ parallel to Mq does not appear 

in this expression. Again, this result does not depend on the initial choice 
of parametrization. One need not choose a parametrization by hand that 
enforces z^uq^ = 0 ; no matter the choice, only the perpendicular part plays 
a role in the field equations. 

The quantity 6T^'^{x;zo) is equal to i’^T^'^(x; zq), where 
is introduced in Appendix |A.2[ and is the Lie derivative introduced in 
Sec.[^ In fact, for any vector T^T^‘'(x; z) is given by Eq. (214) with the 
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replacement —>• and Uq ^ u^. This quantity is useful when considering 
gauge transformations in the self-consistent approximation. Also useful is 
the ordinary Lie derivative of taking similar steps as above, one finds 


C^T^^ix; z) = -mj 


5{x, z) 
(215) 


where = P ^'and .^|| = . The sum of the two Lie derivatives 

yields the simple result 


(£^ -I- £^)T^''{x; z) = -m 


u u u' v 
9^'9u''^ u 



5{x, z)dT'. 


(216) 


A ^8{x, z) term signals that the mass m is displaced from z^ by an 
amount the lack of any Vf^S{x,z) term in Eq. (216) signals that the 
displacements due to the two derivatives cancel one another, leaving the 
mass m moving on 


B.2 Metric perturbation 

In the Lorenz gauge, the first-order self-consistent field, given some global 
boundary conditions, is given by 

hl^{x-z) = 4m j G^ufi'u'U^'u'^'dT', (217) 

where is the Green’s function that comports with the global bound¬ 

ary conditions. I wish to expand this about z^ = Zq to obtain something of 
the form 


^o) + (!^^h]^u{x] zo, zi) -h O(e^). (218) 

There are two methods available for achieving this: directly, following steps 
analogous to those in the previous section; or by making use of the result of 
the previous section. 

Here I adopt the second method. Noting that 6h}^y{x\ zq, zi) = £vhj^j^(x; zq), 
that = £yTl^‘^(x; Zo), and that £y commutes with derivatives acting 

at X, we have 

E;^i,l£vh^](x;zo) = £yP^i,lh^](x; zq) = -167r£yT^^(x; zq). (219) 
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Hence 


= £^h\^{x]ZQ) = 4 / G[x'] zq, zi)dV'. (220) 


From Eq. (214), this evaluates to 


■^0) ^l) - 4?71 / f 2G^j/^/j//'Ug + G y!U q Uq •2]^_|_ j ^Tq. (221) 

•^ 70 ^ ' 


B.2.1 Gauge condition 

It is worth examining how /i^^(x;z) and 6h^i^{x‘, zq, zi) contribute to the 
Lorenz gauge condition. Those contributions are easily found by invoking 
the identity V^G^y^/y' = [9], where is a Green’s function for 

the vector wave equation G\V^ = and both Green’s functions must satisfy 
the same boundary conditions. Performing a trace-reversal on Eq. (217), 


taking the divergence of the result, using the Green’s-function identity, and 
integrating by parts yields 


D 


V^hl,p{x\z) = A / Gaa’jp{rnu^)dT'. 


( 222 ) 


(Earlier results in this section assumed constant m, but here I momentarily 
leave it arbitrary for generality.) The contribution to the gauge condition 
is determined entirely by ^ and the acceleration of z^. In a Gralla-Wald 
expansion, one has zq) = 0, from which one can read off ^ = 0 

and Og =0. In a self-consistent expansion, one instead has V^[e/i)^^(x; z) + 
e^L^/ 3 (x; z)] = O(e^) [or more precisely, Eq. (26)], which determines Eq. ([^. 


Doing the same with Eq. (221) yields 


V^dh^^(x;zo,zi) = 4m / Go 




i± 




+ R^ 


' 7 ^ 5 ' 


dTQ. 


(223) 


The contribution to the gauge condition is determined entirely by the ac¬ 
celeration of the deviation from Zq (together with the geodesic-deviation 
term). In a Gralla-Wald expansion, zq, zi) is included in 20 ), 

and V^L^^(x;zo) = 0 determines Eq. ( [^ . 
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B.3 Local expansion and identification of mass dipole mo¬ 
ment 


In Sec. |4.1t I showed that the mass dipole moment of the object creates a 


term (96) (and contributes to the term (94)) in the object’s skeletal stress- 
energy. By comparing that result to Eq. (214), we can make the identifica¬ 


tion M* = mz\, exactly as concluded elsewhere in the paper. Since the two 
stress-energy tensors are the same, it follows that zq, 2 : 1 ) includes the 

entire contribution to coming from the mass dipole moment. 

Here, I complete the circle by performing a local expansion of Eq. ( |221 ) 
near Zq and showing zq,zi) reproduces the order-1/r^ and 1/r terms 

for the mass dipole seed solution in Fermi normal coordinates. The method 
of local expansion is elaborated in Ref. [9]. In particular, I follow steps 
analogous to those in Sec. 23.2 of that reference. To avoid belabouring the 
details, here I provide only the barest sketch; I leave it to the interested 
reader to fill in the gaps using the tools of Ref. [9] . 

The starting point is the Hadamard decomposition of the retarded Green’s 
function, Here (5+(cr) is a Dirac 

delta function supported on the past light cone of x, and 0+(—cr) is a Heav¬ 
iside step function supported in the interior of that light cone. a{x, x') is 
one-half the squared geodesic distance from x' to x, such that it vanishes 
when a null geodesic connects the two points. 

From this starting point, the final result is obtained by a two-step pro¬ 
cess: (i) evaluating the integrals over Tq by changing the integration variable 
to a, and (ii) writing the retarded distance from x' to x in terms of the Fermi 


radial coordinate r. The results for the two terms in Eq. (221) are 


8m 


where and 




ip.' Dz 


8m 


dr: 


j-dTo = —+ 0(r°), 


, f ^ u' 1/ 'y' 1 I 2 m 

4m / ^0 ^ 


zfriaS^u + C)(r°). 


(224) 


(225) 


The components in Fermi coordinates are 


dhlt = 
dhl = 


’'I'H 


+ O(r0), 


dmii 


r 

2mz\ n. 


+ O(r0), 

+O(r0). 


(226a) 

(226b) 

(226c) 
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Here we see that with the identification M* = mz\, the tt and ab compo¬ 
nents are precisely those in Eq. (105), and the ta component is precisely the 
term proportional to M* in Eq. (106). In other words, the integral (224) 
corresponds to the mass dipole moment’s contribution to the monopole seed 
zq, 6m), while the integral (225) corresponds to the mass-dipole seed 
hl^fix-,zo,M). 


C Identities for gauge transformations of curva¬ 
ture tensors 


Let A\g\ be a tensor of any rank constructed from a metric g. (To stream¬ 
line the presentation, I adopt index-free notation throughout most of this 
appendix.) Now dehne 

1 dA 

6'^A[fi,...,fn] = dAi • • • ^ + • • • + A„/n] (227) 

This tensor is linear in each of its arguments /i,...,/„; it is also sym¬ 
metric in them. In the case that all the arguments are the same, we have 
5'^A[h,... ,h] = ^-^A[g+\h]\^^^, the piece of A[g+h] containing precisely 
n factors of h and its derivatives. 

The following identities are easily proved by writing Lie derivatives as 
ordinary derivatives: 

C^A[g]=5A[C^g\, (228) 

\ClA\g] = \6A[Clg] + 6^A[C^g,C^g], (229) 

C^5A[h] = 6A[C^h] + 25‘^A[C^g, h]. (230) 

Note that 6'^ A[C^g,h] = 6“^A[h, C^g] = ^ (^5"^A[h, C^g] + A[C^g,h]). As an 

example, if A is the Ricci tensor, then 


= (231) 

^C^R^p\g\ = (232) 

£^^6R^ii[}i\ = 6Rp^u\C^h\ -\- 26“^R^p\h, C^g\, (233) 


where I have restored indices to avoid confusion with the Ricci scalar. 

To establish Eq. (228), one can write the metric as a function of a param¬ 
eter A along the flow generated by ^ and then perform a Taylor expansion: 


^iM9\ 


^A 

dX 




A=0 


A=0 


5A 


dAl^=o 


6A[C^g]. (234) 
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Similarly, to establish Eq. (229), one can write 








A=0 




dg 

dX^ 


A=0 


\5A[Clg]+5‘^A[C^g,C^g], 


A=0 


(235a) 

(235b) 


and to establish Eq. ( |230[ ) , one can write g as a function of parameters (A, e), 
where h = ^\ „, and then write 

at le=0’ 

(236a) 

A=e=0 


= JmA 


g{X, 0) + 0) 


dXde 


A 


g(0, 0) + A^(0,0) + e^(0,0) + Ae^(0, 0) 


= 6A[C^h] + 26^A[h,C^g]- 


X=t=0 

(236b) 
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